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NUMERICAL  MODELING  OF  AURORA 
Volume  I of  the  Calendar  Year  1975 
Annual  Report  to  the  Defense  Nuclear  Agency 

Section  1 
INTRODUCTION 

At  the  core  of  the  modern  defensive  systems  are  lines  of 
communication  serving  as  intelligence  links  for  the  initiation  and 
direction  of  countermeasures.  Confidence  in  a communication  system's 
immunity  to  disruption  is  a function  of  one's  understanding  of  the 
large  scale  changes  that  can  be  induced  in  the  communications  environ- 
ment by  both  natural  and  other  means,  such  as  a high  altitude  nuclear 
explosion  HANE). 

A HANE  causes  ionospheric  disturbances  from  prompt  gamma-ray, 
x-ray,  UV,  optical,  or  IR  emissions,  from  late  time  motion  of  ionized 
debris  and  emissions  from  that  debris.  A complex  series  of  physical 
and  chemical  reactions  can  be  initiated  with  consequent  enhanced 
ionization  and  perturbation  of  the  £unbient  electric  and  magnetic 
fields.  The  program  of  the  NRL  Plasma  Dynamics  Branch  is  designed  to 
produce  reliable  models  for  predicting  ionospheric  disturbances  produced 
by  HANEs . 

Confidence  in  a numerical  model  derives  from  the  model's 
ability  to  produce  results  consistent  with  a suitably  complete  and 
accurate  observational  data  base.  However,  the  observational  data  for 
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the  response  of  the  ionosphere  to  nuclear  events  is  insufficient  and 
the  augmentation  of  the  data  is  precluded  by  treaty  restrictions  on 
atmospheric  testing.  Under  these  circinnstances , what  is  needed  is  a 
natural  source  of  large  scale  ionospheric  disturbances  for  which  suffi- 
cient observations  have  been  made.  la  addition  there  must  be  some 
parallels  in  the  physics  of  the  natural  event  and  that  of  a HANE . ilie 
disturbances  responsible  for  the  aurora  meet  these  requirements . 

HANF.s  and  auroras  are  characterized  by  energetic  charged 
particles  streaming  along  geomagnetic  field  lines  and  interacting  with 
the  ambient  atmospheric  constituents.  Particle  energies,  degrees  of 
ionization  and  heating  rates  differ  between  the  two  events,  but  ample 
theoretical  correlations  exist  to  make  one's  capability  for  modeling  the 
aurora  a good  measure  of  one's  capability  for  modeling  the  response  of 
the  ionosphere  to  a HANE . 

Auroral  models  also  help  to  provide  a link  between  the  model- 
ing of  the  late  time  development  of  scintillation  sources  from  HANE 
induced  disturbances  and  scintillation  observations  (e.g.  the  Wideband 
Satellite  Experiment}.  HANE  phenomenology  models  provide  as  output 
the  basic  or  initial  state  of  the  ionospheric  plasma  i.e.  electric 
fields,  currents  and  plasma  densities'  in  the  immediate  aftermath  of  a 
nuclear  event  and  this  in  turn  serves  as  the  input  for  programs  which 
follow  the  development  of  irregularities  in  those  regions  where  the 
plasma  is  unstable.  Tlie  auroral  models  are  analogous  to  the  HANE 
phenomenology  models  in  that  their  output  is  the  basic  state  of  the 
ionosphere,  but  for  an  auroral  disturbance.  The  auroral  precipitation 
patterns  which  serve  as  input  for  auroral  models  can  be  supplied  by 
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s^itellite  observations.  'Ilie  plasma  densities,  currents  and  electric 
fields  which  the  model  computes  are  also  observable  by  satellite  and 
ground  based  radar  and  thus  a check  on  the  auroral  model  is  available. 

An  auroral  model  in  which  confidence  has  been  established  can, 
in  conjunction  with  observations  on  the  ionospheric  basic  state,  serve 
as  a check  on  the  programs  used  to  model  both  HANE  and  aurora-induced 
irregularities.  Observation  of  the  ionospheric  parameters  necessary  as 
input  for  the  progr;ims  which  model  the  formation  of  irregular  regions 
are,  for  practical  reasons,  limited  in  coverage.  Much  scintillation 
data  will  be  taken  in  regions  for  which  no  observations  of  the  basic 
ionospheric  state  have  been  made.  Auroral  models,  in  which  the  output 
quantities  are  computed  self-consistently,  can  be  used  to  determine  the 
missing  data  points,  thereby  extending  the  scope  of  the  scintillation 
observations  . A more  important  consequence  of  this  derives  from  the 
fact  that  the  model  can  provide  the  basic  state  parameters  over  much  of 
the  auroral  oval.  Smaller  scale  computer  programs  then  are  used  to 
follow  the  development  of  irregularities  in  regions  where  plasma 
instabilities  prevail.  These  are  the  same  codes  that  are  used  in  the 
HANE  program.  The  output  of  these  codes  serves  as  input  to  signal 
propagation  codes  which  predict  scintillation  effects.  Comparisons  of 
the  computed  and  observed  scintillations  provides  a measure  of  the 
reliability  of  the  programs  used  to  model  ionospheric  irregularities. 

NRL  has  continued  to  progress  along  three  lines  of  investiga- 
tion which  form  the  core  of  the  auroral  modeling  programs;  (l ) particle 
interactions,  (2)  chemical  reactions,  and  (5)  changes  in  ambient 
ionospheric  currents  and  fields.  The  proton  deposition  code  PRODEP  has 


been  improved  by  the  addition  of  O2  and  O to  the  original  code  which 
had  provided  only  for  the  presence  of  Ng.  The  influence  of  the  geo- 
magnetic field  in  proton  energy  degradation  has  also  been  included. 

Other  work  on  particle  interactions  has  included  a study  of  back- 
scattered  electrons  which  has  confirmed  the  accuracy  of  the  electron 
deposition  code  for  large  angle  scattering.  A summary  of  the  year's 
work  on  particle  interaction  appears  in  section  2 of  this  volume.  The 
work  on  particle  deposition  in  sections  2 and  3 was  done  early  in  CY  75 
and  was  stopped  in  FY  under  the  guidance  of  DNA.  Note  that  this 
annual  report  covers  CY  75* 

Section  h,  on  auroral  chemistry,  contains  a critical  analysis 
of  NO  densities  observed  in  auroral  arcs.  This  study  demonstrates  the 
importance  of  having  an  accurate  knowledge  of  even  the  minor  constitu- 
ents in  the  disturbed  atmosphere.  Finally  we  conclude  with  a 
presentation  of  improvements  in  the  auroral  electrodynamic  models. 

The  importance  of  anomalous  plasma  processes  on  structure  of  electric 
fields  and  current  systems  has  been  studied  in  detail,  as  the  struc- 
ture of  the  fields  and  currents  are  important  to  ionospheric  plasma 
stability  and  scintillation  phenomena.  Also,  we  present  improved 
model  results  for  the  midnight  region  of  the  auroral  oval.  It  is 
these  types  of  results  which  will  be  of  great  use  in  interpreting  Wide- 
band Satellite  Data. 

The  work  performed  at  NRL  during  the  past  year  is  summarized 
in  the  remainder  of  this  volume.  Much  of  this  work  has  been  previously 
reported  in  NRL  Memo  Reports,  at  several  symposia,  in  technical  jour- 
nals, and  at  DNA  sponsored  meetings.  The  principal  contributors  in  each 
technical  area  are  listed  as  authors  of  the  section  describing  that  work. 


Section  2 


IMPROVED  MODEL  TO  STUDY  PROTON 
DEPOSITION  IN  THE  ATMOSPHERE 

J . E . Rogerson  and  J . Davis 

The  deposition  of  energetic  particles  into  the  earth's 
atmosphere  causes  an  abundance  of  interesting  phenomena.  In  addition 
to  the  spectacular  visual  display  produced  by  auroral  emissions,  the 
auroras  provide  a unique  means  for  studying  the  complex  of  mechanisms 
occurring  in  the  atmosphere.  Knowledge  of  the  causes  of  these  effects 
is  important  in  assessing  the  influence  the  disturbed  atmosphere  has  on 
communication  systems  and  their  design.  Also,  since  enhanced  ioniza- 
tion is  generated  in  the  polar  region,  changes  in  the  atmosphere's 
characteristics  can  be  investigated  during  PGA  events. 

Man"  idles  of  the  auroral  electron  deposition^"®  and 
proton  dep  has  been  reported.  In  an  earlier  report,  a 

simpliPi  based  on  a continuous-slowlng-down  approximation 

(CSDA)  was  described.  This  model  contained  only  one  atmospheric  con- 
stituent, Ng,  and  neglected  the  earth's  magnetic  field.  In  this 
revised  model,  the  number  of  atmospheric  constituents  is  extended  by 
including  N^,  0^,  and  0.  By  assuming  a dipole  model  of  the  geomagnetic 
field,  an  expression  for  the  differential  element  of  arc  length  trav- 

1 0 "7 

ersed  by  a charged  particle  spiraling  around  a magnetic  field  line 
is  incorporated  into  the  CSDA  to  produce  a more  realistic  proton  energy 
degradation.  The  charge  exchange  formalism  as  given  by  Allison^®  is 
generalized  to  a three  component  atmosphere  and  a charge-equilibrated 
flux  of  protons  and  hydrogen  atoms  is  generated  from  which  atmospheric 
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excitation  and  ionization  processes  can  be  calculated.  Secondary 
electron  effects  are  treated  as  before^^.  Finally^  total  ionization 
results  are  presented  and  compared  with  results  of  the  previous  model 
and  with  similar  calculations  by  Father  and  Burrows®,  who  used  a differ- 
ent technique. 


2.1  THE  MODEL 


A charged  particle  moving  in  a magnetic  field  follows  a 
spiral  path  about  the  field  lines  due  to  the  Lorentz  force.  The  arc 
length  of  this  spiral  path  depends  on  the  particle's  charge,  velocity, 
mass,  and  pitch  angle  with  respect  to  the  given  field  line.  An  equation 
for  the  differential  element  of  arc  length  of  the  spiral  path  can  be 
given  as 


ds 


" t (i-x; 


■<  -s  in^  0!  '4“5x 


dr 
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where  ds  = increment  of  arc  length, 

dr  = increment  of  geocentric  distance, 

r = equatorial  radius  of  the  magnetic  field  line  fL  shell), 

^ I 

X = r/r^, 

= equatorial  pitch  angle. 

Eq  . f2.l)  assumes  a dipole  magnetic  field  given  by 
B (r,x)  = M v^i+j  sin^X  / 
where  M = dipole  moment 

\ = magnetic  latitude. 

>,  and  r are  related  by 

r = r cos^  > . ',2.5' 

e 


(2.^4) 


The  first  adiabatic  invariant 
sin^  (y/B  = constant 
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is  also  included,  where  a is  the  pitch  angle. 

Kq . , .1  contains  a singularity  at  x = 1 , or  r = r^ . This 

singularity  corresponds  to  a point  in  the  equatorial  plane.  For 
auroral  studies,  lilgh  latitude  regions  are  of  primary  interest,  and  it 
is  therefore  reasonable  to  avoid  this  point.  A singularity  also  occurs 
for 

.X  - sin‘"’  (V  h-~.x  = 0 , 

which  corresponds  to  a mirror  point  ’u  = ’~/2).  This  second  singularity 
corresponds  to  an  initially  dowmward-moving  particle  being  reflected 
upward  due  to  the  converging  magnetic  field  lines.  Such  a particle 
does  not  penetrate  as  deeply  as  it  might  if  it  had  a smaller  pitch 
angle;  hence,  it  does  not  produce  as  much  ionization  as  it  otherwise 
might.  Also,  the  treatment  of  "mirroring''  particles  introduces  fur- 
ther complexities  into  the  calculation.  For  this  reason,  and  because 
the  aim  is  to  stud-'  the  maximum  perturbations  created  in  the  atmosphere, 
this  singularity  is  also  avoided. 

By  assuming  a spherical  earth  of  radius  a,  Eq . '2-1  can  be 
taken  as  a function  of  altitude  z by  setting  r = a+z ; then 
X = (a+z)  / r^ . For  a given  altitude  difference,  z,  " ^2  ' 
can  be  numerically  Integrated  from  z^  to  z^  to  give  the  arc  length 
traveled  by  the  particle  in  going  from  to  z^ . Then,  by  using  the 
CSDA  as  described  in  Reference  15,  the  proton  energy  degradation  is 
computed  as  follows; 

AE  = Nj.  (z^)  Lp  (E,z^)  AS  , {2.b) 

where  N^(z^)  is  the  total  atmospheric  density  at  z^,  Lp(E,z^)  is  the 
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total  proton  loss  function  at  for  energy  E,  and  AE^  and  AS  is  the 
arc  length  traversed  in  going  from  z^  to  z^ . 

liie  individual  loss  function  for  protons  in  0^,  and  0 
are  adapted  from  analytic  expressions  given  by  Green  and  Peterson^'^. 

At  altitude  z,  the  total  loss  function  is  computed  from 

N . ^z  ) 

-E  sV  Si '2-«) 

i=l 

where  i is  the  species  index. 

At  a given  altitude  and  energy,  a generalized  charge  equili- 
brated flux  is  calculated  based  on  Allison's^®  charge  exchange 
formalism.  The  generalized  equilibrated  fraction  of  protons  and 
hydrogen  atoms  are  given  by 


S n'-  'z)  Oq^’-I'E) 

^ N^(z)  [aQ^'-(E)  + a^Q^'E  )] 
i 

(2.7) 


XI  z)  a^Q^;E) 

F^(E,z)  = ^ ^ 

= 1 - F^!'E,z'  , 

where  the  index  i runs  over  the  three  species  N^,  0^,  and  0.  is 

the  cross  section  for  charge  exchange  of  a proton  with  a neutral 
'creating  a hydrogen  atom),  and  is  the  cross  section  for  electron 

stripping  of  a hydrogen  atom  by  a neutral  (creating  a proton  again). 
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For  a given  process  a involving  proton  impact  on  species  i,  a volume 


production  rate  at  altitude  z for  a given  proton  flux  <t>  (E ) is 
Ppi  /Fi(E,z)  (E)  ^ 'E)  dE, 


(2.8) 


and  similarly  for  hydrogen  atom  impact  processes  with  F^^E,z)  replaced 
by  F ;E  ,z  ) . 

To  describe  the  creation  of  secondary  electrons  of  kinetic 
energy  T by  direct  ionization  processes,  an  energy-differential  cross 
section  of  the  form 

S(E,T)  = da(E)/dT 


a.  (E) 

ion 


r tan 


ri  + rvr  ] 


-1 


'2.9) 


as  given  by  Edgar^  Miles,  and  Green^'^  is  used.  For  E in  keV,  and 


T in  eV, 
max  ^ 


m 

^ * 
^ 1 

n 


T = 2.18  E. 
max 


This  is  equivalent  to  a form  given  by  Banks,  Chappell,  and  Nagy^  for 
secondary  electron  production  from  electron  impact  ionization.  By 
using  the  P values  given  in  Reference  3 and  the  value  of  22  eV  given 
by  Reference  for  proton  impact  ionization  of  N^,  r values  for  proton 
impact  ionization  of  0^  and  0 can  be  obtained  by  scaling.  The  same 
S(E,T)  was  used  for  both  proton  and  hydrogen  atom  impact  on  the  same 
target  gas  because  of  insufficient  data.  Secondary  electrons  created 
by  the  charge  exchange  cycle  are  treated  exactly  as  in  Reference  l4; 
i .e . , 
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E,T  = 


CVC  It* 


F.  at 


"c/1' 


T 

■<>  ni.TX 


. 1 


^^?here  a cycle 
each  species . 
electrons  are 
Julienne® . 


- -■  + - . anti  I't  is  taken  to  be  ■ eV  for 

Excitation  and  ionization  by  secondary  and  higlier  order 
treated  by  the  metliod  described  by  Dalgarno  and  used  by 


2.2  CROSS  SECTIONS 

The  principal  sources  for  H and  impact  cross  section  data 
are  the  review  articles  by  McNeal  and  Birely”  j'  The  charge  exchange 
cross  section  for  0^  is  taken  from  Stebbings  , et  al  for  lower  ener- 
gies and  by  doubling  the  0 charge  exchange  cross  section  as  given  oy 
Allison^®  for  the  higher  energies.  The  electron  stripping  reaction  of 
H by  Og  is  obtained  from  doubling  the  0 cross-section  data  given  by 
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Allison  for  the  same  process.  The  charge  exchange  cross  section  for 
0 was  taken  from  Neynaber,  et  al."^®  for  energies  below  one  kev  and  froia 
McNeal  and  Birely^^  for  higher  energies.  Many  of  these  cross  sections 
do  not  extend  below  two  to  four  keV  in  energy  and  were  extrapolated  to 
lower  energies . 

The  input  cross-section  data  points  for  proton,  hydrogen 
atom,  and  electron  impact  processes  are  fitted  to  the  points  on  a chosen 
energy  grid  E by  a cubic  spline  representation  as  a function  of  In  E. 

In  this  way,  only  a few  experimental  points  are  needed  for  a good 
representation,  and  analytic  fits  to  cross  sections  are  not  required. 

The  cross  sections  are  extrapolated  to  high  energies  by  using  the 
appropriate  asymptotic  forms. 


'ITie  electron  cross  sections  are  summarized  in  Table  2.1 


2.3  RESlfLTS 

In  this  section,  some  results  of  the  revised  proton  deposi- 
tion model  are  presented  and  compared  with  the  previous  model^^  and 
with  similar  calculations  by  Eather  and  Burrows  (hereafter  denoted  by 
EB  ). 

Figure  2.1  shows  a comparison  of  the  energy  degradation  of  an 
incident  100  kev  proton  in  the  two  NRL  proton  deposition  models  for  a 
pure  atmosphere;  in  these  calculations,  the  I963  Jacchia^**  model 

Cl 

atmosphere  is  used.  Results  for  pitch  angles  of  zero  degrees  and 
for  about  45  degrees  are  showii  when  the  earth's  magnetic  field  is 
taken  into  account.  It  can  be  clearly  seen  that  the  effect  of  the 
magnetic  field  is  to  greatly  increase  energy  degradation  and  thus 
reduce  the  penetration  depth  into  the  atmosphere. 

Figure  2.2  shows  the  total  number  of  N^"*”  ions  produced  as  a 
function  of  altitude  for  these  same  100  keV  proton  deposition  calcula- 
tions. At  higher  altitudes  there  is  very  little  difference,  but  at 
lower  altitudes,  fewer  ions  are  created  when  the  earth's  magnetic  field 
is  included  because  of  the  reduced  penetration  into  the  denser  regions. 

Figures  2. 3"^. 6 show  comparisons  of  some  results  of  the 
revised  NRL  model  with  similar  calculations  by  EB . In  these  computa- 
tions, the  1961  Chamberlin^^  model  atmosphere  is  used,  and  N^,  0^,  and 
0 are  included  as  constituents.  In  the  EB  results,  monoenergetic 
protons  with  isotropic  pitch  angle  distributions  are  considered,  whereas 
monoenergetic  protons  with  a single  pitch  angle  are  treated  in  the  NRL 


results.  Depositions  in  the  polar  regions  are  considered  in  both 
cases.  Figures  2.J-2.6  present  comparisons  of  the  total  ionization 
produced  per  unit  volume  versus  altitude  at  several  incident  energies. 

It  can  be  seen  in  these  Figures  that  the  EB  results  at  the 
lower  energies  can  be  matched  very  well  at  the  lower  altitudes  by  vary- 
ing the  pitch  angle.  At  higher  altitudes^  there  is  some  disagreement, 
especially  for  the  1000  keV  calculation.  The  large  differences  for  a 
1000  keV  proton  are  suggestive  of  a significant  difference  in  the  models 
used  by  EB  and  by  NRL.  Further  study  of  References  9 i*nd  25  indicated 
that  the  cross  section  for  the  process 

h'*'  + X + x'*'  + e,  '2.12^ 

i.e.,  direct  ionization  by  protons,  was  not  included  in  the  EB  calcula- 
tions. At  lower  energies,  charge  exchange  is  the  dominant  process  for 
ionization  by  protons,  and  the  equilibrium  fraction  of  hydrogen  atoms 
is  much  larger  than  for  protons;  hence,  this  is  a reasonable  approxi- 
mation at  lower  energies.  At  higher  energies  charge  exchange  cross 
sections  are  greatly  reduced,  the  equilibrium  fraction  of  protons  is 
large,  and  the  contribution  from  direct  ionization  is  very  significant. 

To  check  whether  this  was  indeed  the  case,  the  depositions 
were  reconputed  with  the  cross  sections  for  the  processes  given  by 
Eq . (2.12/  set  equal  to  zero.  By  then  redefining  the  total  number  of 
charged  particles  at  a given  altitude  as  the  sum  of  the  ^2^^ 
o"*"  ions  plus  the  number  of  protons  generated  from  t}ie  charge  exchange 
cycle,  a new  set  of  results  was  obtained.  These  results  are  sho\vai  in 
F igures  2 .7-2  .10  . 
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It  is  seen  that  £or  energies  of  100  keV  or  less,  the  agree- 
ment with  EB  is  still  good.  In  Figure  2-7,  adding  the  charge  exchange 


protons  greatly  overestimates  the  peak  of  the  curve  by  about  50  pet 
cent,  whereas  agreement  is  very  good  without  this  contribution.  As 
noted  in  the  previous  section,  however,  the  cross  sections  were 
extrapolated  to  very  low  energies,  and  this  is  probably  the  major 
source  of  the  discrepancy.  (The  results  in  Figure  2.7  for  total  ions 
only  arc  the  same  as  the  results  in  Figure  2-5.  This  follows  from  the 
fact  that  at  5 keV,  the  equilibrium  fraction  of  hydrogen  atoms  is 
approximately  unity;  i.e.,  direct  ionization  by  protons  is  effectively 
turned  off . ) <\greement  is  good  for  either  calculation  at  JO  keV  and 
100  keV , as  might  be  expected  since  this  is  the  energy  regime  where 
charge  exchange  is  dominant.  The  notably  better  agreement  occurs  at 
1000  keV,  as  may  be  seen  by  comparing  Figures  2.6  and  2.10.  While  the 
agreement  in  Figure  2.10  is  not  totally  satisfactory,  it  is  much  better 
than  that  represented  by  Figure  2.6  and  appears  to  verify  the  assumption 
that  the  proton  direct  ionization  process  is  omitted  by  EB . The 
remaining  discrepancies  are  probably  due  to  differences  in  cross- 
section  data. 

The  results  presented  thus  far  have  involved  only  monoener- 
getic  proton  depositions.  The  program  also  treats  a distribution  of 
proton  energies.  .'\n  input  flux  with  the  energy  distribution  shown  in 
Figure  2.''.  1 was  inserted  into  the  program,  and  the  results  for  total 
ionization  produced  as  a function  of  altitude  are  shown  in  Figure  2.12. 
In  this  calculation,  the  Chamberlin  I96I  atmosphere  was  again  used, 
although  it  is  possible  to  use  any  given  model  atmosphere. 
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Results  given  in  this  report  have  dealt  only  with  ionization 
produced  by  proton  depositions.  However,  it  is  possible  to  calculate 
the  creation  of  other  excited  states  or  emission  lines  by  putting  in 
the  relevant  cross  sections.  For  example,  the  calculation  of  the 
Ng  + angstrom  emission  line  can  be  calculated  by  using  the 

McNeal-Birely^  ' data  for  this  cross  section. 

The  deposition  model  can  be  applied  to  deposition  of  energetic 
ions  in  debris  particles  such  as  that  associated  with  STARFISH.  In  the 
coming  year  work  will  be  done  in  this  regard. 
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Fig.  2.1  — Comparison  of  100  keV  proton  energy  degradation 
atmosphere  with  and  without  the  geomagnetic  field  Bg.  a 
pitch  angle  with  respect  to  the  field  lines. 
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Fig.  2.2  — Number  of  ions  produced  vs  altitude  with  and  without  the 
geomagnetic  field  Bg  by  deposition  of  a 100  keV  proton  in  a pure  S2 
atmosphere,  a is  the  pitch  angle. 
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P’ig.  2.5  — Number  of  ions  produced  by  a 100  keV  proton  as  a function  of 
altitude  in  a Chamberlin  1961  model  atmosphere 
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FiS.  2.7  — Recalculation  of  the  ionization  produced  by  a 3 keV  proton  in  a 
1961  Chamberlin  model  atmosphere  when  the  cross  section  Oijj  for  direct 
ionization  by  protons  is  set  equal  to  zero 
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2.8  Recalculation  of  ihe  ionization  produced  by  a 30  keV  proton  in 
a 1901  rhamix'rlin  model  atmosphere  when  the  cross  section  Oj,,  for  direct 
ionization  by  protons  is  set  equal  to  zero 
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Fig.  2.10  — Recalculation  of  the  ionization  produced  by  a 1000  keV  proton 
in  a 1961  Chamberlin  model  atmosphere  when  the  cross  section  Ojp  for  direct 
ionization  by  protons  is  set  equal  to  zero 
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Section  ;j 

TRANSPORT  TECHNIQUES  FOR  DESCRIBING 
THE  ELASTIC  SCATTERING  OF  ENERGETIC  ELECTRONS 

D.  J.  Strickland  and  I.  B.  Bernstein 
In  the  past  year,  a code  has  been  developed  at  the  Naval 
Research  Laboratory  which  solves  an  equation  of  transfer  for  energetic 
electrons  The  code  is  being  applied  to  auroral  studies  and  deter- 
mines, for  a given  incident  auroral  electron  spectrum,  the  electron 
flux  as  a function  of  pitch  angle,  energy,  and  altitude.  The  back- 
scattered  flux  near  l8o°  is  very  sensitive  to  the  detailed  angular 
dependence  of  the  cross  sections.  In  spite  of  the  fact  that  the  cross 
sections  for  large  scattering  angles  are  orders  of  magnitude  smaller 
than  for  angles  close  to  zero  degrees,  changes  in  the  cross  sections  at 
large  scattering  angles  produce  noticeable  changes  in  backscatter.  To 
investigate  such  angular  behavior  in  detail,  we  have  applied  the  program 
to  a simpler  problem,  that  of  purely  elastic  scattering  by  one  constit- 
uent. The  accuracy  of  the  results  is  easier  to  determine  for  this 
case . 

We  begin  by  discussing  the  equation  of  transfer  or  Boltzmann 
equation  for  elastic  scattering.  This  is  followed  by  introducing  a 
cutoff  angle  for  a single  scattering  in  order  to  study  the  importance 
of  large  angle  scattering.  Next,  the  Fokker-Planck  equation  for  elas- 
tic scattering  is  introduced  since  it  is  derivable  from  the  equation  of 
transfer.  Both  numerical  and  analytic  solutions  have  been  obtained  for 
the  equation  of  transfer,  and  analytic  solutions  have  been  obtained  for 
the  Fokker-Planck  equation.  A subsection  is  devoted  to  discussing  these 


methods  of  solution.  Finally,  selected  results  for  the  two  equations 
are  presented  followed  by  our  conclusions. 

3 .1  METHODS 

We  wish  to  describe  electron  distribution  functions  wtiich  may 
possess  strong  anisotropies.  'The  Boltzmann  equation  provides  a starting 
point  for  such  an  exercise  but  in  its  most  general  form  is  not  practical 
for  most  applications  due  to  its  complexity.  There  are  several  simpli- 
fied versions  of  this  equation  which  have  had  extensive  use  over  tlu- 
years,  especially  in  the  area  of  plasma  physics.  Of  these  we  shall 
consider  two  here,  the  Boltzmann  equation  for  infinitely  massive 
scatterers  and  the  counterpart  Fokker-Planck  equation,  both  for  thi 
case  of  elastic  scattering  alone. 

We  begin  with  the  Boltzmann  equation  governing  the  distribu- 
tion function  f for  elastic  scattering  by  infinitely  heavy  scatterers: 

,t  ^ 

where  n is  the  density  of  scatterers,  a is  the  differential  cross 
section,  and  is  the  scattering  angle  see  Fig.  *.l  . Consider  a 
steady  state  in  the  absence  of  electromagnetic  fields  where  f depends 
on  the  space  coordinate  z alone.  Eq . '^.1  then  reduces  to 

U,  ^ vf(z,v)  = dn  o . v,'^  )v'‘ f z ,y' ) - f z ,v  ' 

where  u,  is  the  cosine  of  the  pitch  angle  u = v-e^/v  . Define  the 
total  cross  section 
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Fig.  3.1  — Geometrj'  for  the  scatU'ring  of  a particle  from  direction 
n'  to  n referenced  to  fixed  direction  z 
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and  the  differential  scattering  depth 


dr  = n-i  dz  . 
o 


Assimiing  azimuthal  symmetry  in  f,  Eq . (5.2)  then  yields 

U ~ du'  R(p,,  u'  ' 

which  is  conventionally  termed  the  equation  of  transfer. 
The  kernel  R is 
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A<t>  ■j(v^'=')/c 
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Consider  the  case  of  the  screened  Coulomb  cross  section 
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2 4 
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(l-COS  '=!  + 21]  J‘ 
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with  screening  parameter  T]  . 

On  introducing 


cos  = v-v'/v^  = u,u'  + '^1-u-^  cos 


R becomes 


p-  M 2Ti:i  + h)(l  + 2T-uu'\  ^ 

’ '■(u-u'F  + 4Tia-uM.')  + 


(5-7) 


(5-8) 


(5-9) 


Consider  next  a cutoff  angle  ^ such  that  for  larger  values 
of  the  scattering  angle  ,0  is  set  to  zero.  The  parameter  enables 

us  to  examine  the  dependence  of  f on  large  angle  scattering.  In  terms 
of  this  modified  cross  section,  R is 
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^ ^ tan 

a + b cos^  V a‘  -b‘ 

max 


Va^-b^’  max 
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where  a = 1 + 2T]  - u,n,%  b = /l-u"  = cos  and  is 

the  maximum  allowed  azimuthal  angle  for  given  U,  and  |j,'.  See 
Appendix  A for  further  details;. 

To  solve  Eq . (5*5  the  range  in  u,  has  been  divided  into  N 
zones,  and  the  flux  has  been  given  a piece-wise  linear  dependence  over 
the  grid.  The  matrix  representation  of  Eq . (5 *5)  is  then 
df. 

u.  ^ = FA.  .f.  . '^.11 ; 

t IJ  J 

Examples  of  the  dependence  of  R on  U and  u,^  are  given  in 

Fig.  3.2  for  a value  of  T|  of  10"“*.  Such  a value  would  be  appropriate 
for  electrons  in  the  100  kilovolt  range  scattering  off  a light  nucleus. 
The  sharply-peaked  behavior  of  R at  small  scattering  angles  suggests 
that  one  make  a Fokker-Planck  expansion  of  the  collision  integral  in 
Eq.  (3.2)  (see,  e .g . Bethe  et  al.^®  and  Rosenbluth  et  al."^).  The 
resulting  equation  [see  Appendix  B]  is 


^ ^ ^ (1  - 

2a^  Ap,  |_  Ap 


(3.12) 
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= d'^  sin  ^ (1  - cos  '=')a(v,'^) 
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is  the  momentum  transfer  cross  section.  For  the  screened  coulomb 
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In  our  finite  difference  approximation  to  Eq . (5  •12),  the 
right  hand  side  is  represented  by 


g^Aq^  [(l-q|.l)(F._^- 
0 ^ 

■f.  ) - (l-q^  l)(f.-f._^  )1 

1 ^1+  i 1 1+1  J 

(5-15) 

at  interior  points,  and  by 

CT  Aq 

fq  = 1) 

(3-16) 
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^ All  N-1 

0 

1 — 1 
tl 

:± 

(3-17) 

where  the  domain  in  u.  has  been  divided  into  N equal  intervals  and 
u.  is  the  value  of  u,  half-way  between  q,.  and  q..  We  have  further 
multiplied  the  terms  corresponding  to  q = ± 1 by  - Aq/4)  for 
reasons  explained  at  the  end  of  this  section. 

Both  numerical  and  analytic  methods  have  been  used  to  solve 
the  matrix  equations.  Specifically,  both  methods  have  been  applied  to 
the  equation  of  transfer,  and  the  analytic  technique  to  the  Fokker- 
Planck  equation. 

Physically,  for  a well  posed  problem  one  must  stipulate  the 
distribution  function  for  those  particles  incident  on  the  system; 
collisions  then  determine  the  properties  of  those  particles  that  are 
emerging.  For  the  case  here  considered  this  corresponds  to  stipulating 
f o,q)  for  1 2 q > 0 and  f(r^,q)  for  0 > q S -I-  In  the  work  here 
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presented  the  choice  for  all  cases  was  f(Tg,u)  = 0 , 0 > p,  & - 1> 
corresponding  to  an  absorbing  boundary. 

The  solutions  obtained  by  numerical  integration  are  based  on 
an  iteration  scheme  as  follows;  First^  for  0 < s 

to  zero  and  using  a second  order  predictor-corrector  scheme,  Eq . '''..11 
is  integrated  forward  to  t = the  lower  boundary)  for  0 < ^ 1 

and  in  the  process  the  results  are  stored  on  a pre-specified  grid  in 
-.  Eq . j .11 ) is  then  integrated  backwards  to  t = 0 for  0 ''  p S -1 
using  the  values  of  f^(r)  just  computed  for  < p^  < 1,  and  just  as 
with  the  forward  Integrations,  the  results  are  stored  on  the  same  grid 
in  T.  Fon^ard  and  reverse  integrations  proceed  until  the  values  of 
f^(r)  do  not  change  appreciably.  In  practice,  a technique  is  used 
which  projects  the  iterated  values  of  f^  to  their  converged  values  in 
only  three  or  four  iterations . 

We  have  used  a standard  eigenvalue  method  to  obtain  the 
analytic  solutions  of  the  difference  differential  equations.  Assuming 
f^  can  be  represented  by  'ii^e  we  solve  the  characteristic  equation 


T. 

j 


(i  = 1,N) 


(3.18) 


leading  to 


fi(T)  = E C^lTi^e'^'/  (f  = 1,N) 


(3-19) 


where  \ is  the  1*"^  eigenvalue,  * . is  the  eigenfunction  at  u , 

^ 1 0 


and  C is  a coefficient  determined  by  boundary  conditions, 

I 
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By  inspection,  one  can  verity  that  the  following  is  also  a 
solution  to  either  the  equation  of  transfer  or  the  Fokker-Planck 
equation; 


. f . 
1 


O'  + 3 


20 


where  a and  9 are  constants  for  further  details  see  Bethe  et  al.  . 

It  turns  out  that  two  of  the  eigensolutions  to  the  discre- 
tized equations  serve  to  represent  this  solution,  namely  those  with 
eigenvalue  which  go  to  zero  as  N ->  oo.  In  order  to  minimize  errors  for 

T » 1 associated  with  approximating  the  exact  solution  (Eq . 5 -20 ) by 
o 

a sum  of  exponentials,  it  is  convenient  to  suppress  the  two  eigen- 
solutions corresponding  to  the  smallest  eigenvalues  and  write 


: . = a + 9 p,  - ^ )+  J: 
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If  one  seeks  eigensolutions  (p,  )e  to  Eq . 3-12),  they 
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from  wnich  it  follows  tnat 
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45 

Hence,  eigensolutions  of  the  differential  equation  carry  no  current. 

To  enhance  the  correspondence  between  eigensolutions  to  the 
Fokker-Planck  matrix  equation  and  Eq . (5.I2),  the  former  are  made  to 
satisfy  condition  (5.23)  in  the  form 


.5-* 
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(J.2^ ; 


N-l 

y,  + 2 r Vi  + yN  = ^ 

i = 2 

wlicre  y^  within  the  sum  is  given  by  our  finite  difference  approxima- 
tion to  the  right  hand  side  of  Eq . (5.22)  ''see  Eq  . (5.15))-  Equation 
(5.24)  then  specifies  y^  and  y^^  and  leads  to  the  replacement  of  Eqs  . 

' .l6)  and  (5.17  ) by 

-^)(Vi  - <n) - -1' 

5.2  RESULTS 

In  this  section,  we  present  solutions  of  the  equation  of 
transfer  and  of  the  Fokker-Planck  equation  for  a variety  of  conditions. 
Values  of  10”^  and  10"“^  for  the  screening  parameter  are  considered, 
which  would  apply  to  kilovolt  electrons  scattering  off  light  nuclei. 

A wide  range  of  total  scattering  depths  is  assumed,  from  values  of 
5 to  1000.  Values  of  the  cutoff  angle  of  50°  and  110°  have 

been  used  in  the  equation  of  transfer.  Finally,  incident  fluxes  with 
both  isotropic  and  Gaussian  dependences  have  been  assumed. 

Figs.  5.5  and  5.4  demonstrate  the  sensitivity  of  the  back- 
scattered  flux  to  ct(P',71)  at  large  pitch  angles.  The  solutions  were 
obtained  with  the  equation  of  transfer.  The  screening  parameter  is 
1 for  results  in  Figs.  5.5  and  10"°  for  those  in  Fig.  5.4,  thus 
providing  for  much  stronger  forward  scattering  for  the  first  set  of 
results.  The  total  scattering  depth  is  100  for  both  systems.  In 
each  figure,  results  are  given  at  t = 0 and  50  for  values  of  180  , 
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PITCH  ANGLE 


Fig.  3.3  — Solutions  to  the  equation  of  transfer  at  the  upper  boundar>’  and 
halfway  through  the  medium  for  cutoff  angles  <r)^  of  180°,  30°,  and  15° 
and  for  t}  = lO""* 
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30°  and  15° • We  remind  the  reader  that  scattering  is  restricted  to 


values  less  than  or  equal  to  the  chosen  value  of  . The  backscattered 
flux  near  l8o°  does  decrease  significantly  as  large  angle  scattering 
is  restricted  in  spite  of  the  fact  that  the  values  of  the  cross  sec- 
tions being  affected  are  extremely  small  (see  Fig.  3*2  for  ’ = 11”'^). 
The  results  can  be  explained  in  terms  of  the  effective  turning  point 
for  the  particle.  Once  the  particle  does  turn,  its  probability  of 
directly  escaping  through  the  upper  boundary  goes  as  exp  --  where 
T is  the  depth  from  this  turning  point  to  the  boundary.  A small  change 
in  T arising  from  varying  produces  a much  larger  change  in  the 
escape  probability,  which  then  reflects  itself  in  the  backscattered 
flux.  The  results  in  Figs.  3 >5  ^nd  .k  lead  to  an  important  point  -- 
namely  that  a detailed  knowledge  of  the  scattering  cross  section  is 
essential  to  accurately  predict  backscatter  for  systems  with  not 
much  greater  than  1.  The  Fokker-Planck  solutions  exhibit  a similar 
sensitivity  to  through  changes  that  occur  in  the  momentum  transfer 
cross  section. 

Figs.  3-5  arid  3 -6  give  fluxes  at  t = 0 for  T]  = 10"^  and  for 

5 ^ T ^ 1000.  Since  all  results  refer  to  the  same  value  of  T],  the 
o 

ratio  of  the  total  number  of  scatterers  for  any  two  situations  is  sim- 
ply the  ratio  of  the  corresponding  values  of  . The  incident  flux  is 
isotropic  in  Fig.  5 "5  and  assumes  a Gaussian  dependence 
'"exp  (-1  "u,  )^/ .25  ) J in  Fig.  3.6.  For  each  case,  solutions  have  been 
obtained  to  both  the  integral  equation  and  its  counterpart  Fokker- 
Planck  equation  and  are  presented  in  these  figures.  As  x^  increases. 


5- -3 


4Kim‘ 


the  albedo  or  ratio  of  total  backscatter  to  total  incident  flux 
approaches  unity.  Results  for  both  equations  are  in  good  agreement 
for  systems  with  t >■  100.  For  such  cases,  enough  scatterers  are 
present  for  multiple  small  angle  scattering  to  be  effective.  With 
decreasing  , however,  single  large  angle  scattering  becomes 
relatively  more  important  and  is  not  adequately  described  by  the 
Fokker-Planck  equation,  which  underestimates  the  true  amount  of  back- 
scatter . 

The  variation  of  the  flux  with  - is  shown  in  Fig.  5.7  for 
T|  = 1 = 10,  and  isotropic  incidence.  Again,  solutions  to  both 

equations  are  given.  Since  we  are  considering  a rather  thin  system, 
the  incident  flux  experiences  very  little  attenuation.  The  solution 
at  " = 1 illustrates  the  lower  boundary  condition:  zero  flux  from 

below.  A second  example  of  the  variation  of  the  flux  with  - is  shown 
in  Fig.  5 -8  • Here,  the  system  is  ten  times  thicker  and  the  incident 
flux  has  a Gaussian  dependence  on  p, . Fokker-Planck  solutions  are  not 
included  because  of  the  number  of  curves  already  presented.  Their 
behavior  can  be  fairly  well  estimated  from  considering  the  previous 
figures.  The  imposed  minimum  at  t = 0 due  to  the  incident  profile 
rapidly  disappears  with  increasing  -.  For  t > 5>  the  flux  becomes  a 
monotonically  decreasing  function  with  p. 


X 

Fig.  3.9  — Relationship  between  incident  velocity  v^and  resultant  velocity 
v'  = V + A.  The  case  treated  is  for  elastic  scattering,  i.e.  v'  = v. 
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Table  3.1.(Q/2a^!  from  Eq . (3.1^+)  versus  Screening  Parameter  T] 
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Section  4 

AURORAL  CHEMISTRY  --NO 

E.  Hyman,  D.  J.  Strickland, 

P.  S.  Julienne  and  D.  F.  Strobel 

A most  intriguing  aeronomical  problem  in  recent  years  is  the 
anomalously  large  enhancements  of  NO  concentration  that  have  teen 
inferred  from  measured  rN0'*”]/[02'*']  ratios^^""^®;  [Swider  and  Narcisi, 
1970;  1974]  and  directly  measured®®  in  auroral  arcs.  'The  analysis  of 
Jones  and  Rees®^  suggests  that  a measurement  of  the  '"N0"'']/[02''"]  ratio 
is  not  necessarily  an  accurate  means  to  infer  the  NO  density  at  the 
time  of  the  measurement.)  The  available  evidence  suggests  that  these 
enhancements  are  spatially  and  temporally  restricted  to  auroral  arcs. 
There  is  considerable  variability  in  the  peak  NO  densities;  for  example, 
Swider  and  Narcisi'^®  infer  [NO]  ~ 6x10®  cm"®  at  110  km  whereas  Zipf 
et  al.®°  measured  3 cm"®  at  120  km.  It  should  be  noted  that 
other  measurements  on  the  latter  rocket  flight  suggest  ''NO]  was  only 
1-5x10^°  cm"®  at  110  km  as  compared  with  the  mass  spectrometer  result 
of  2x10^^  cm"®.®^ 

Attempts  to  deduce  NO  concentrations  from  studies  of  the 
deactivation  of  N2(A®E^)  molecules®^  and  from  NO  y band  emission®® 
are  not  conclusive  since  the  production  mechanisms  are  not  properly 
understood.  In  the  former  case  the  experimental  cross  sections  adopted 
in  the  study  were  too  high.®"^  In  the  latter  case  it  has  not  been 
demonstrated  that  auroral  y band  emission  is  indicative  of  ambient  NO 
concentrations  . 
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The  Zipf  et  al.  measurements  are  thus  the  most  perplexing 
observations  to  explain  and  appear  to  violate  the  principle  of  energy 
consert'ation  if  precipitating  auroral  primaries  are  the  major  energy 
source.’’^  Previous  theoretical  studies  of  auroral  events  have 
addressed  a broad  range  of  auroral  processes  but  have  not  specifically 
tackled  the  formidable  problem  of  auroral  NO.  In  this  paper  we  con- 
fine our  attention  to  auroral  NO  and  investigate  the  applicability  and 
limitations  of  conventional  midlatitude  ionospheric  chemistry  and 
certain  speculative  chemical  schemes  to  address  the  question  of  how 
much  NO  can  be  produced  during  an  auroral  event  and  how  chemically 
this  production  can  be  maximized.  We  conclude  that  while  localized 
enhancements  of  '"NO']  to  ~ 10^°  cm”^  are  possible  in  a narrow  altitude 
range,  special  and  improbable  circumstances  are  required;  for  example, 
vibrationally  hot  0^ . 


4.1  CONVENTIONAL  CHEMISTRY 


First  we  discuss  auroral  NO  production  on  the  basis  of  con- 
ventional midlatitude  chemistry,  where  odd  nitrogen  is  formed  by 
dissociation  and  ionization  of  N^  with  9*8  arid  15*6  eV  required, 
respectively.  In  the  case  of  N^^  production  subsequent  ion-molecule 
reactions  are  needed  to  break  the  N^  bond, but  these  reactions  are  very 
exothermic.  In  either  event  considerable  energy  (>  lOeV  ' must  be 
expended  . 


In  the  altitude  range  of  interest,  100-120  km,  N^  is  the 
dominant  constituent  and  consequently  the  precipitating  energetic 
electrons  will  dissipate  most  of  their  energy  by  dissociating  and 
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ionizing  N . There  are  three  principal  channels  that  must  be  con- 
sidered and  are  illustrated  in  Figure  4.1.  A direct  ionization  channel 
results  in  the  production  of  which  either  reacts  with  0 to  form 

NO^  and  N ' D or  with  0 to  recycle  • For  typical  fluxes  and 
auroral  composition,  li of  the  particle  flux  energy  is  expended  to 
produce  odd  nitrogen  by  this  channel.  If  FNO]  ~ [0]  then  N,  ^ would 
preferentially  react  with  NO  to  give  NO^  + and  cut  off  odd  nitrogen 
production . 


also  dissociates  in  the  manner  described  by  Oran  et  al . 
to  produce  N'^D)  + N('^S).  Approximately  26^^  of  the  total  energy  is 
consumed  in  this  channel.  In  addition  11*^  of  the  particle  energy  is 
used  to  dissociativelv  ionize  N and  form  N^  and  N.  The  electronic 

C 

state  of  the  N atom  is  not  kno^^n  but  presumed,  without  serious  conse- 

4 + 

quences,  to  be  in  the  ground  S state.  The  N will  react  with  0^  to 
produce  either  NO"*^  or  N(^D).  According  to  studies  by  Rusch  et  al.^® 
and  Strobel  et  al.^^,  recent  odd  nitrogen  data  requi'-es  that  NO"*^ 
dissociatively  recombine  to  form  principally  N(^D).  Thus,  the  princi- 
pal source  of  N(^S  is  the  pure  dissociation  channel. 

The  remaining  ~ 50*!?  of  the  energy  is  used  to  ionize  0 and 

-j- 

Og,  to  produce  radiation,  etc.  Subsequent  reactions  of  0^  and  0 do 
not  produce  odd  nitrogen,  except  the  minor  reaction  path 
o'*’  + Ng  NO"^  + N('^S).  In  the  region  100-120  km  o'*'  reacts  preferen- 
tially with  0^  unless  either  0^  is  translationally  hot  and/or  is 

vibrationally  excited.  The  measured  0^  concentrations  preclude  the 
*4" 

rapid  removal  of  0 by  a fast  reaction  path  and  imply  that  0 has  less 
than  0.4  eV  of  translational  energy  and  that  the  v s 2 vibrational 


levels  of  N have  a Boltzmann  distribution  at  the  kinetic  temperature. 


In  summary,  the  precipitating  energetic  electrons  expend 
~ of  their  energy  in  the  production  of  odd  nitrogen.  It  takes 

I 

about  10  eV  to  dissociate  and  about  2(  eV  to  produce  a pair  of  odd 
nitrogens,  witich  is  more  efficient  than  the  55  eV  required  to  produce 
an  ion  pair.  Thus,  while  odd  nitrogen  is  efficiently  produced,  sub- 
stantial energy  (lOeV  - is  required  to  initiate  its  production.  The 
chemistry  outlined  in  Fig.  4.1  results  in  the  production  of  1.9  N “D 
atoms  for  each  N ' atom.  The  molecular  proce.sses  involved  in  N atom 
production  are  not  sensitive  to  the  environmental  conditions  in  the 
auroral  zone  in  a manner  that  would  alter  the  N(^D)  to  N ''S)  production 
ratio. 

For  N(^D)  and  N(‘*S)  there  are  two  major  reaction  paths:  a 

fast  reaction  with  NO  and  a slower  reaction  with  0^ . In  the  case  of 
N(^D)  the  reaction  paths  are  comparable  when  [N0]/[0  1 ~ 0.1.  Thus, 
until  NO  is  built  up  to  large  concentrations  N(^D  mostly  reacts  with 
0^  to  produce  more  NO.  For  N('^S  > it  is  well  known  that  its  reaction 
with  0^  involves  a substantial  activation  energy  (~  6 k cal  and  at 
room  temperature  is  very  slow  ~ 10"'®  cm^  sec”^.  The  reaction  paths 
of  N(^S)  are  comparable  when  [NOj/fO^]  ~ 4x10"®;  i.e.,  for  NO  concen- 
trations be^ow  ambient  values.  Large  buildups  of  NO  are  precluded  in 
a short  period  of  time  under  normal  circumstances  as  much  of  the  NO 
produced  by  N(^D)  + 0^  will  be  destroyed  by  N('^S)  + NO . It  is  obvious 
that  to  maximize  the  auroral  NO  production  rate  the  appropriate  chemis- 
try must  favor  conversion  of  NC^S)  into  NO  rather  than  N^ . Accordingly 
the  reaction  NC’s)  + 0.  must  proceed  at  a rapid  rate.  The  small  steric 
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factor  for  this  reaction  eliminates  the  possibility  of  hot  N''^S)  atoms 
reacting  with  0 before  being  themalized  by  collisions  with  the  more 

abundant  N molecules.  Electronic  excitation  of  0^  to  state  does 

not  help  since  the  reaction  N(‘’s)  + 0 (^A  ) -»  NO  + 0 is  slow 

cm  ^ sec”  , Clark  and  Wayne'®"'.  An  examination  of  the  reac- 

tion trajectory  suggests  that  0 vibrational  energy  should  be  available 
to  surmount  the  energy  barrier.  For  sufficiently  hot  0^  v N '"'S  could 
react  at  a rate  approaching  2x10“^^  cm ^ sec”^  and  ’"N0^/'’0^]  ~ .15 

would  be  the  condition  for  equal  probability  of  reaction  with  0^  and 
NO  . 

Even  with  vibrationally  hot  0^  in  an  auroral  arc,  the  buildup 
of  NO  would  terminate  when  '"N0]/'’02]  ??  0.1,  since  any  additional  N ''s 
and  N(^D)  atom  production  would  destroy  rather  than  increase  the  NO 
concentration.  At  120  km  we  estimate  that  NO  could  not  exceed 
5x10^  cm  ® in  contradiction  to  the  mass  spectrometer  measurements  of 
Zipf  et  al.’^  This  upper  limit  TnO]  estimate  would  have  to  be  lowered 
if  large  amounts  of  0^  were  dissociated  as  the  spectrometer  results  of 
Zipf  et  al.  suggest.  It  should  be  emphasized  that  given  sufficient 
time  10^  seconds  for  a typical  auroral  arc)  the  NO  concentration 

will  build  up  to  0.1^  if  the  N '^D)  production  rate  exceeds  that  of 
N(''s)  and  if  transport  effects  are  not  included. 

k.2  THE  MODEL 

To  make  these  estimates  quantitative  we  have  constructed  an 
idealized  model  to  describe  the  buildup  of  NO  during  auroral  events. 

The  model  is  described  in  detail  in  a technical  report  available  from 
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the  authors.  Production  rates  of  ions,  excited  states,  N 'S  , 

N ■ I)  , optical  emission,  etc.,  are  j^enerated  by  an  energetic  electron 
deposition  code  described  in  a technical  report  by  Strickland  et  al.  ' 
The  severity  of  the  auroral  NO  problem  allows  us  to  neglect  transr-ort 
processes.  Our  model  results  may  yield  larger  NO  densities  than 
measured  as  a consequence  of  this  assumption.  For  initial  values  of 
odd  nitrogen  densities,  the  model  of  Strobel®"  was  adopted;  the 
results  are  insensitive  to  this  choice.  The  time  dependent  continuity 
equations  are  integrated  for  all  important  constituents  vhose  chemical 
time  constant  exceeds  one  second.  Chemical  equilibrium  is  assumed  for 
constituents  with  chemical  time  constants  less  than  one  second.  The 
important  rate  constants  are  tabulated  in  Table  h .1 . 

The  model  was  applied  to  describe  the  remarkable  auroral  arc 
probed  by  Zipf  et  al.  The  energy  deposition  rate,  illustrated  in 
Fig.  4*2,  was  chosen  to  reproduce  the  measured  .!  intensity  as  a 

function  of  altitude.  Under  the  most  favorable  conditions,  T ' > 

2 eV,  the  NO  buildup  as  a function  of  time  and  altitude  is  illustrated 
in  Fig.  4.2.  According  to  Donahue^^  this  auroral  arc  existed  previous 
to  the  measurement  for  ~ 4x10^  seconds  . Our  model  predicts  a maximum 
'"NOl  ~ SxlO®  cm"^  with  the  enhancement  localized  in  the  95"H5 
region.  To  reach  the  NO  cutoff  density  ~ O.lFO^]  requires  buildup 
times  » 1-'^  seconds. 

The  sensitivity  of  the  model  to  various  input  parameters  is 
illustrated  in  Fig.  4. 5.  The  NO  buildup  is  most  sensitive  to  the  0^ 
vibrational  temperature  and  the  N(^d)  branching  ratio  for  NO  + e. 
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Fig.  4.2  — Altitude  profiles  of  the  [NOj  buildup  as  a function  of  time 
for  vibrationally  hot  O2,  Ty(02>  = 0.2  eV.  The  energy  deposition  rate 
is  also  illustrated  and  produces  26  kR  of  3914  A radiation  at  100  km. 
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With  '>0<^  N ‘*S  ) production  from  NO^  + e,  the  NO  buildup  is  severely 
limited  by  the  large  chemical  loss  from  the  reaction  N'^S  + NO . How- 
ever, such  a large  N('*S)  production  is  not  consistent  with  midlatitude 
data.^"’®^  Our  results  are  not  sensitive  to  the  selection  of  initial 
NO  densities.  The  sensitivity  of  the  NO  buildup  to  a decreased  0 
content  of  the  model  atmosphere  is  due  to  a reduction  of  odd  nitrogen 
production,  as  more  N^^  reacts  with  0^  rather  than  with  0.  This  reduc- 
tion in  NO  production  exceeds  an  opposite  tendency  due  to  reduced 
N(^D'  quenching  by  0.  An  increase  in  the  0 content  produces  less 
pronounced  '"NO"’  changes. 

4.5  CRITIQUE  OF  MODEL 


The  model  can  describe  a rapid  buildup  of  NO  if  the  0^ 
vibrational  temperature,  T^(O^),  is  large,  but  the  sources  and  sinks 
of  0^  vibrational  energy  must  be  compatible  with  high  T^(O^)  in 
auroral  arcs.  Measurements  at  l600-3500°  K indicate  that  0 atoms  are 
efficient  quenchers  of  0^  vibrational  energy.^  Extrapolating  their 


results  down  to  300°  K,  we  obtain  a quenching  rate  of  ~ 10 


-IS 


sec"*  for  Og (v=l ) . Two  possible  sources  of  0^  vibrational  excitation 
are  the  precipitating  electrons  and  collisions  of  0^  with  translation- 
ally  hot  ions.  In  the  former  case  less  than  5^  of  the  energy 
deposition  is  channeled  into  0^  vibrational  energy.  Trans lat ional ly 
hot  ions,  in  particular  0^^,  may  be  able  in  collisions  to  efficiently 
excite  02(v'.  This  latter  process  would  be  physically  similar  to 
Joule  heating  and  maximum  excitation  would  occur  at  ~ 1?5  km  with 
negligible  excitation  below  llC  km  where  the  maximum  ’’NO"'  enhancement 
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occurs  ^cf.  reference  52  and  Fig.  h .2~' ■ ’fhe  incorrect  altitude  varia- 
tion and  the  observed  ion  composition  preclude  translationally  hot  ions 
as  the  source  of  0^(v).  li/hereas  the  precipitating  electrons  could 
excite  0^  in  the  appropriate  altitude  region,  its  source  strength  is 
not  potent  enough  to  overcome  the  severe  quenching  by  0 atoms  to  give 
high  T^(0^''.  High  T^(O^)  is  possible  only  if  the  0 quenching  rate  were 
comparable  to  the  slow  0^  and  quenching  rate  at  -jPO' ' K.  We  do  note 

that  Reid  and  Withbroe®^  suggest  0^  vibrational  enhancement  in  the 
thermosphere  from  observations  of  solar  UV  radiation.  A room  tempera- 
tuie  measurement  of  the  0 quenching  rate  would  be  decisive.  On  the 
basis  of  our  present  knowledge,  it  is  unlikely  that  the  0 quenching 
rate  is  ~ 10”^®  cm®  sec“~  at  500°  X and  hence  high  T^'O^)  is  improbable 
in  auroral  arcs  where  large  '‘NO*  enhancements  are  inferred. 


SPECULATIVE  CHEMISTRY 

In  light  of  our  previous  discussion  our  search  for  potent 
auroral  NO  sources  must  become  more  desperate.  In  conventional 
chemistry  approximately  10  eV  of  energy  must  be  expended  to  initiate 
odd  nitrogen  production.  The  reaction  ^2  ^ *^2  ^ only  endo- 

thermic by  ~ 2 eV  and  a potential  factor  of  10  Increase  in  the  NO 
production  efficiency  is  possible.  There  is  no  experimental  evidence 
to  indicate  that  this  reaction  occurs  even  with  vibrational ly  excited 


Og  and/or  N^ • 


But  if  it  were  possible  to  form  NO  by  this  reaction. 


there  would  be  no  chemical  cutoff  in  the  production  of  NO  as  present 


in  the  conventional  chemical  cycle. 
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The  reiictlon  + 0 ->  NO  + N could  also  be  considered  with 


race  constant  ~ 10  exp 


2Xp  ^ 


-■:)3ooo\  3 


T / 


cm"^  sec"  . To  compete  with  conven- 


tional  NO  sources,  K ~ 5x10"^'*  sec"'  or  T ~ 2000'  K.  If  N^  were 
vibrationally  hot  and  this  energy  were  available  to  overcome  the  acti- 
vation  barrier,  0 atoms  would  quench  (v J with  a rate  ~ ^xlO  cm 
sec"^,®'^  much  faster  than  the  rate  at  which  N and  NO  would  be  produced. 
The  calculations  of  Schunk  and  Hays®^  support  the  conclusion  that  N^ 
must  have  a Boltzmann  distribution  for  vibrational  energy  in  the  E- 
region  during  auroral  bombardment.  Anharmonicity  effects  at  higher 
vibrational  levels  would  not  negate  their  conclusions.  Reaction 
N^  + 0 cannot  be  the  source  of  anomalous  auroral  TNO]  . 


4.5  DRIZZLE  REGION 

The  available  evidence  suggests  that  the  auroral  zone  under- 
goes continual  bombardment  by  soft  energetic  electrons  and  only 
occasionally  do  spectacular  events  occur.®®  The  continual  bombardment 
of  the  auroral  regions  will  lead  to  an  enhanced  background  NO  concen- 
tration. With  the  model,  an  attempt  was  made  to  estimate  this 
concentration  in  the  absence  of  transport  for  an  electron  deposition 
rate  that  produces  1 kR  of  the  auroral  green  line  (5577  ^ ),  its 
threshold  for  detection.  This  deposition  rate  is  equal  to  0.05  of  that 
in  Fig.  4.2.  A continual  increase  in  the  NO  concentration  with  time  is 
obtained  until  the  NO  cutoff  density  is  reached  at  a time  » 10® 
seconds.  After  10®  seconds,  vertical  transport  effects  will  become 
important  and  severely  limit  the  further  NO  buildup.  However,  peak  NO 
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concentration  ~ 10‘‘  cm“'^  would  be  possible  if  the  above  electron 
deposition  rate  were  constant  witli  time.  Witli  such  concentrations, 
the  TNO]  inferred  by  Donahue  et  al."^^  and  Swider  and  Narcisi'"'®  from 
ion  composition  measurements  would  be  possible  without  exotic 
chemistry,  although  it  could  not  explain  the  peculiar  kink  in  their 
deduced  [NO]  profile.  In  the  case  of  the  mass  spectrometer  measure- 
ments of  large  ’’NO]  this  source  of  NO  fails  by  at  least  one  order  of 
magnitude . 

i;  .6  CONCLUDING  REMARKS 

On  the  basis  of  our  analysis  we  conclude  that  ''NO''  enhance- 
ments in  excess  of  0.l[02]  are  exceedingly  improbable  in  auroral  arcs. 
As  a consequence,  we  are  very  skeptical  of  mass  spectrometer  measure- 
ments that  exceed  this  limit.  Any  TNOI  enhancement  past  this  limit 
would  be  inconsistent  with  the  chemistry  discussed  here  with  one 
exception,  NO  formation  by  N^  • We  emphasize  that  there  is  no 

experimental  support  for  its  occurrence.  All  other  chemical  schemes 
have  the  same  difficulty;  the  NO  buildup  stops  abruptly  when  [NO] 
reaches  ~ O.lfO^],  because  NO  is  not  directly  formed  from  N^.  The  N^ 
must  be  first  dissociated;  then,  the  N atoms  are  oxidized.  A maximum 
NO  buildup  rate  is  achieved  under  these  conditions  only  if  the  reac- 
tion N('’^S)  + 0^  ->  no  + 0 is  fast.  The  most  plausible  mechanism  is 
vibrationally  hot  0^  but  limited  data  at  shock  tube  temperatures 
suggest  that  0 quenching  is  too  severe  to  permit  a departure  from 
thermal  equilibrium.  A room  temperature  measurement  of  O^iv  = l)  + 

0 0^(o)  + 0 is  desirable. 
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NO  BUILDUP  AT  100  KM 
(O2 ) IN  ev 


;[N0] 

amb  ~ 2.0[N0],„b 


'l0]-0.5[0] 


0,0166 


+e'  -*0(^P)  + 0.5[n('’s)  + N(2d)1 


TIME  (sec) 

Fig.  4.3  — The  [NO]  buildup  for  various  parameters  as  a function  of 
time  at  100  km  where  peak  deposition  occurs.  For  the  standard  run, 
denoted  by  0.0166,  only  N(2D)  is  formed  from  NO'*’  + e,  Ty(02)  = 
0.0166  eV  = kinetic  temperature,  and  [O]  = 7 X 10*®  cm~®  at  120 
km.  The  label  of  the  curve  indicates  the  parameter  varied.  The  labels, 
0.2  and  0.6,  denote  the  values  of  enhanced  T^,(0.,). 
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Section  5 

EFFECTS  OF  ANOMALOUS  RESISTIVITY  ON 
AURORAL  BIRKEL/\ND  CURRENT  SYSTEMS 

Joel  A.  Redder 

In  recent  years  it  has  become  widely  recognized  that  geomag- 
netic field-aligned,  Birkeland,  electric  currents  are  of  considerable 

importance  to  the  physics  of  the  high-latitude  ionosphere  and  magneto-  ' 

i 

sphere.  Although  their  existence  in  association  with  aurora  was  j 

I 

orginallv  proposed  by  Birkeland,  their  importance  only  became  li 

1 

evident  as  a result  of  the  great  number  of  measurements  of  magnetic  | 

I 

I 

field  perturbations  in  the  auroral  regions  by  rockets  and  satellites. 

The  rocket  measurements*^'  ^ have  shown  clearly  the  association  of 
Birkeland  currents  with  auroral  arcs  and  bands.  Whereas,  the  satellite 

measurements have  shown  that  the  field-aligned  currents  occur  over  ' 

a much  wider  region  associated  within  the  whole  auroral  oval. 

The  Birkeland  currents  in  the  auroral  region  arise  from  two 
different  mechanisms.  Energetic  plasma  pressure  gradients  iu  the 
magnetosphere  can  lead  to  field-aligned  currents  as  discussed  by 
Vasyliunas  Additionally,  the  ionosphere  can  drive  field-aligned 

currents  caused  by  an  imposed  electric  field  in  the  presence  of  con- 
ductivity gradients  or  by  an  ionospheric  conductivity  in  the  presence 

of  electric  field  gradients.  These  Birkeland  currents  couple  the  ■ 

magnetospheric  and  ionospheric  electric  fields  and  to  some  extent  must 

Vs 

modify  the  spatial  distribution  and  magnitude  of  each. 

Models  for  the  field-aligned  currents  and  electric  fields  i 

have  been  developed  for  auroral  arcs  and  for  the  auroral  region  as  a 
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whole.  Bostron/''  presented  a model  for  electric  fields  and  current  in 
arcs  making  two  limiting  assumptions  concerning  the  magnitude  of  the 
Birkeland  current.  In  one  case  the  Birkeland  current  was  allowed  to 
flow  freely  while  in  the  other,  it  was  taken  to  be  negligible. 

Atkinson®'  also  presented  a model  for  an  auroral  arc  in  which  the 
return  ionosphere  electron  current  was  strongly  limited.  In  this 
paper  results  are  presented  which  are  similar  to  theirs.  Models  for 
the  electric  fields  and  currents  in  the  auroral  regions  and  lower 
latitudes  as  well  have  been  developed  by  Wolf,®  Jaggi  and  Wolf,®* 

Iwasaki  and  Nishida,®^  Karlson,®®  and  others.  In  most  of  this  work 
the  conductivity  of  the  plasma  along  the  geomagnetic  field  has  been 
taken  to  be  large,  and  therefore  they  show  no  effects  of  field-aligned 
resistance  on  the  derived  Birkeland  currents  and  electric  fields.  One 
exception  to  this  is  the  model  of  Coroniti  and  Kennel®"^  for  the  polari- 
zation of  the  westward  electrojet  in  the  midnight  auroral  oval,  caused 
by  anomalous  resistance  to  field-aligned  current.  However,  the  idea  of 
anomalous  resistance  to  field-aligned  currents  did  not  arise  as  a 
result  of  comparing  model  results  to  data  for  auroral  electric  fields 
and  current,  but  rather,  from  a desire  to  explain  the  acceleration  of 
auroral  energetic  particles. 

Swift®®  proposed  the  existence  of  the  ion  acoustic  insta- 
bility in  the  auroral  plasma  which  could  lead  to  anomalous 
resistivity,  a field-aligned  electric  potential,  and  the  acceleration 
of  auroral  electrons.  Subsequently,  Kindel  and  Kennel®®  showed  that 
for  auroral  plasma  parameters,  the  ion  cyclotron  instability  had  a 
lower  threshold  and  therefore  would  be  a more  likely  mechanism  for  the 
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anomalous  resistivity.  Furthermore,  the  threshold  for  this  current 
driven  instability  is  well  within  the  range  of  Birkeland  current 
densities  inferred  from  measurements.  A further  study  of  the  ion 
cyclotron  instability  by  Palmadesso  et  al.®^  demonstrates  the  self- 
quenching character  of  the  instability  which  results  in  marginal 
instability  and  considerably  lower  resistivity.  In  addition  to  the 
current  driven  plasma  instabilities  mentioned  above  Papadopoulos  and 
Coffey®®^®®  have  proposed  a beam  driven  parametric  instability  which 
also  may  provide  resistance  to  auroral  Birkeland  currents.  In  this 
mechanism,  a precipitating  electron  beam  drives  high-frequency  plasma 
waves  which  interact  parametrically  to  drive  lower  frequency  waves  and 
plasma  density  fluctuations.  This  parametric  Instability  has  the 
virtue  that  in  addition  to  providing  anomalous  resistivity,  it  also 
stabilizes  the  beam  plasma  interaction  and  enhances  the  tail  density 
in  the  thermal  electron  distribution  to  agree  with  auroral  particle 

90 

measurements  . 

Mechanisms  other  than  plasma  instabilities  and  anomalous 
resistance  have  been  suggested  which  could  lead  to  field-aligned 
electric  potentials.  It  would  be  expected  that  these  mechanisms  also 
would  greatly  affect  Birkeland  current  systems . One  of  these  is  the 
possible  formation  of  electrostatic  double  layers®^^®®  which  could  lead 
to  a large  potential  drop  over  a small  distance  along  the  magnetic 
field.  The  other  is  the  possible  formation  of  electrostatic  shock 
waves  between  the  magnetosphere  and  ionosphere.  Kan®®  has  discussed 
normal  electrostatic  shocks  while  Swift®"^  has  discussed  the  possible 
occurrence  of  oblique  electrostatic  shocks  and  their  effects  on 
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particles  and  electric  fields. 

In  this  section  we  attempt  to  discuss  in  considerable  detail 
the  effects  of  anomalous  resistance  on  auroral  Birkeland  currents.  We 
will  illustrate  and  delineate  these  effects  by  showing  a rather  large 
selection  of  results  from  numerical  model  studies.  In  the  following 
section  the  important  features  of  the  auroral  model  are  presented.  We 
also  will  present  a simple  model  for  the  ionosphere-magnetosphere 
interaction  with  anomalous  resistivity  to  field-aligned  currents. 

The  section  on  results  will  demonstrate  the  effects  of 
anomalous  resistivity  on  Birkeland  currents.  Results  are  presented  for 
ooth  the  minimum  and  maximum  magnitude  of  resistivity  which  can  reason- 
ably be  expected.  Also  to  be  presented  are  model  results  for  the  two 
most  likely  resistive  mechanisms.  We  show  the  effect  of  treating  the 
precipitating  current  separately  from  the  ionospheric  electron 
Birkeland  current.  These  results  are  discussed  and  some  conclusions 
will  be  drawn. 

Since  the  simple  model  used  to  show  the  effects  of  anomalous 
resistivity  is  probably  not  very  realistic  for  comparison  to  real 
auroral  arcs,  in  the  final  section  we  will  present  and  discuss  the 
results  of  one  model  for  the  midnight  region  of  the  auroral  oval.  The 
complex  current  systems  in  this  model  will  be  discussed  in  some  detail. 
We  will  also  demonstrate  that  the  electric  fields  derived  from  this 
model  are  in  agreement  with  satellite  electric  field  measurements. 
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MATHEMATICAL  MODEL 


To  describe  the  auroral  ionosphere,  we  solve  the  plasma 
continuity  equation 

^ + V • 'nv  ''  = P - on^  , 

At  i ~L 

the  plasma  momentum  equation 

^ 'E  + — V X B > = vv  , 
m c ~i  - X " 

and  the  current  continuity  equation 

V • J = 0 . 

In  Eq . (5.1'  we  have  neglected  vertical  transport  of  the  plasma;  in 
Eq . (5.2),  the  inertial  terms;  and  in  Eq . charging  effects.  In 

Eq.  (5.1),  n,  V,  P,  and  a are  respectively  the  plasma  density,  the  ion 
velocity,  the  plasma  volume  production  rate,  and  the  volume  dissocia- 
tive recombination  rate;  and  the  symbol  ; applies  to  the  geomagnetic 
field.  The  symbols  in  Eq . (5.2),  q,  m,  E,  c,  B,  and  v are  respec- 
tively the  electron  or  ion  charge,  and  electric  field  in  the 
ionosphere,  the  speed  of  light,  the  geomagnetic  field  and  the 
electron-  or  ion-neutral  collision  frequency.  The  current  density, 
in  Eq . (5.5),  '-S  defined  by 

1 = en  1 V - v . ) , ' 5 

where  e is  the  electron  charge  and  v and  v.  are  the  electron  and  ion 

— e —1 

velocities . 

These  equations  are  solved  on  a cartesian  coordinate  system 
with  the  x-axis  directed  north-south;  the  y-axis,  east-west;  and  the 
z-axis,  vertical.  In  this  paper,  for  simplicity,  it  is  assumed  that 
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all  ionospheric  properties  are  constant  along  the  y-dimension.  Addi- 
tionally, the  geomagnetic  field  is  taken  to  be  directed  along  negative 
z.  Although  these  assumptions  are  not  necessary  to  achieve  a valid 
solution,  they  do  make  the  results  much  easier  to  interpret. 

In  order  to  complete  the  set  of  Eqs . (5 -I  - 'y -h  , it  is 
necessary  to  prescribe  the  electric  field.  We  define  E in  terms  of 
the  electric  potential, 

E = , (5-5) 

which  is  defined  by  the  magnetosphere-ionosphere  interaction.  Ohms 
Law; 

+ J||R  . (5-6) 

^ is  the  ionospheric  electric  potential,  the  prescribed  magneto- 
spheric  potential  on  the  same  field  line.  J||  is  the  field-aligned 
current  leaving  or  entering  the  top  of  the  horizontally  conducting 
ionosphere  km  altitude),  and  R is  the  resistivity  integrated 

along  a field  line  from  200  km  altitude  to  the  equatorial  plane  in 
the  magnetosphere.  Eq . (5.6)  implies  that  the  magnetosphere  is  an 
intrinsic  source  of  electric  potential  which  drives  the  ionosphere  and 


Birkeland  current  systems.  Although  taking  the  auroral  magnetosphere 
as  a .urce  of  either  electric  field  or  current  is  not  usual  (see 
discussion  by  Wolf  ),  other  workers  have  made  such  an  assumption.''  ’ 
'■qs  . . - (5.6)  are  a complete  set  and  may  be  solved 

self-consistently  to  yield  the  plasma  density,  electric  field,  and 
rrent  svstem  provided  that  P (in  Eq.  (5-l))  and  and  R (in 
Eq.  . are  specified.  The  production  rate,  P,  was  chosen  to 
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closely  resemble  the  ' emission  data  of  Romick  and  lielon. 

Using  tvpical  values  for  the  recombination  rate  and  the  collision 
frequencies,  the  maximum  Pedersen  conductivity  in  the  auroral  arc  is 
enhanced  a factor  of  .2  over  that  outside,  and  the  [lall  to  Pedersen 
conductivity  ratio  is  2.4.  We  judge  that  these  values  are  tvpical  of 
manv  auroral  arcs,  is  chosen  to  yield  an  electric  field  perpendicu- 

lar to  the  arc,  in  the  x direction,  which  has  a strength  of  2S 
millivolts  per  meter.  ITiis  is  the  simplest  situation  for  the  electric 
field  about  an  auroral  arc,  and  therefore  allows  us  to  more  clearly 
show  the  effects  of  the  resistance  on  the  Birkeland  current.  The 
resistance  function  has  been  chosen  differently  for  each  result  and 
will  be  described  with  the  results. 

To  obtain  a self-consistent  result,  it  is  necessary  to  obtain 
a simultaneous  solution  of  (5-1^  - (5-6).  The  results  presented  here 
were  obtained  as  follows: 

a.  Using  assumed  , obtain  ^ from  (5.6)  and  (5 .5) 
with  (t  = 0 ) = 0 , 

b.  Obtain  v from  '^.2), 

c.  Advance  n in  time  with  (5.l)> 

d.  Obtain  J||  using  (5-5  ^^d  (5.4), 

e.  Compute  E from  (5.6)  and  (5-5), 

f.  Repeat  b - e until  a steady  state  is  achieved. 

The  equations  are  solved  for  the  altitude  range,  95  “ 200  km,  where  the 
perpendicular  ionospheric  currents  occur.  The  results  include  n and  2 
throughout  this  region,  and  and  J||  at  200  km  altitude. 


5-7 


'ajp'*'*; 


RESULTS 


[n  this  section  the  effects  of  anomalous  resistivity  on 
auroral  Birkelancl  currents  are  explored  by  consideration  of  a group 
of  examples.  We  first  consider  the  effect  of  the  magnitude  of  the 
resistivity  on  the  amplitude  of  the  Birkeland  current  density  and  on 
the  electric  field  perturbation.  We  next  consider  the  results  of 
anomalous  resistivity  for  two  of  the  proposed  plasma  mechanisms;  the 
current  driven  ion  cyclotron  instability  and  the  beam  driven  parametric 
instability.  Finally,  results  are  shown  for  a number  of  cases  in  which 
the  precipitating  current  is  resistance  free  and  only  the  return  cur- 
rents experience  resistivity. 

Figure  ‘p.l  shows  the  result  for  the  minimum  amount  of 
resisti-ity  that  can  be  expected  along  an  auroral  field  line.  The 
resistance  was  chosen  at  a constant  value  of  3.5  x 10"^  e.s.u.  This 
value  is  the  Spitzer®'"^  resistance  of  a 1 cm'’  flux  tube  above  200  km  at 
auroral  latitudes.  The  Birkeland  currents  peak  strongly  at  the  edges 
of  the  arc  while  the  electric  field  perturbation  is  negligible.  The 
maximum  current  density  has  the  value  5 x 10^°  el  cm“®  sec'^ 

(8.65  X 10"^  amp  m"^ ) which  is  considerably  larger  than  current 
densities  normally  observed.  This  result  is  quite  similar  to  one 
originally  discussed  by  Bostrom.^® 

Figure  5.2  shows  a contrasting  result  to  Figure  5 • In  this 
case  the  resistance  is  again  constant  but  now  is  the  maximum  value  that 
can  reasonably  be  expected  along  auroral  field  lines.  Its  value  is 
2 X 10^'  larger  than  the  Spitzer  value  and  is  comparable  to  that  used  in 
previous  work  by  Coroniti  and  Kennel.®”^  For  such  a large  value  of 
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ELECTRIC  FIELD 


Fig.  5.1  — Electric  field  and  Birkeland  current  density  for  a very  small 
resistance  to  field-aligned  current.  The  resistance,  R = 0.00035  esu, 
(the  Spitzer  resistance  of  a 1 cm^  flux-tulie  above  200  km  altitude)  is 
spatially  constant.  For  this  small  resistance  the  electric  field  pertur- 
liation  is  small  in  the  auroral  arc  while  the  Birkeland  current  density 
is  very  large,  up  to  8.65  X 10~®  amp 
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ELECTRIC  FIELD 


BIRKELAND  CURRENT 


P’ig.  5.2  — Electric  field  and  Birkeland  current  density  for  a very  large 
resistance.  The  resistance,  7 esu,  is  spatially  constant  and  a factor  of 
2 X lO-*  larger  than  the  Spitzer  value.  For  this  large  resistance  the 
Birkeland  current  density  is  500  times  smaller  than  in  the  previous  case. 
Whereas,  the  electric  field  is  greatly  reduced  inside  the  auroral  arc. 


p 

. ' resistance,  the  results  show  that  the  i5irkeland  current  is  small 

1 .1  X 10^  el  cm"'  sec"^  1 .'('6  x I-"  amp  m”‘  ) and  that  the  electric 
field  is  greatly  reduced  in  the  auroral  arc.  The  Birkeland  current 
density  is  smaller  than  what  is  normally  observed  whereas  the  electric 
field  strength  is  anti-correlated  with  the  arc  as  has  been  observed , 

As  in  the  previous  case,  this  type  of  situation  has  also  been  discussed 

by  Bostrom.^®  I 

The  results  in  Figure  5-3  ate  intermediate  between  those 
shown  in  the  first  two  figures.  In  this  case  the  resistance  has  the  j 

constant  value  0.3  esu.  The  electric  field  and  Birkeland  current  i 

patterms  are  quite  similar  to  those  of  the  previous  case.  But  here,  | 

the  current  density  has  a maximum  of  h x 10"  amp  m""^,  which  is  more  j 

I j 

typical  of  measured  values  . In  all  of  these  cases  where  a constant  | 

. resistance  was  chosen,  the  Birkeland  currents  maximize  at  the  strong 

i I 

density  gradients  on  this  side  of  the  auroral  arc.  On  the  other  hand,  j 

the  electric  field  shows  a minimum  in  the  center  of  the  arc  and  tends  ^ 

to  have  maxima  outside  the  arc  structure. 

For  the  next  result,  we  attempt  to  model  a marginally  unstable 
current  driven  resistance  mechanism.  Marginal  instability  is  a reason- 
able condition  as  the  current  driven  instabilities  are  self-saturating. 

The  marginal  instability  model  assumes  that  the  resistance  on  each 

i 

flux  tube  is  incrementally  increased  above  the  Spitzer  value  at  each 
time  step  when  the  Birkeland  current  density  exceeds  a threshold  of 

3 X 10®  el  cm"''  sec  The  result  is  shown  in  Figure  The  i 

electric  field  is  quite  similar  to  that  in  the  previous  case,  however, 

the  Birkeland  current  is  considerably  different.  The  Birkeland  current  ’ 
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ELECTRIC  FIELD 


Fig.  5.3  — Electric  field  and  Birkeland  current  density  for  an  intermediate 
value  of  resistance.  The  resistance  is  spatially  uniform  and  has  a viUue  of 
0.28  esu.  The  Birkeland  current  density  has  a reasonable  magnitude,  about 
4 X 10~®  amp  m~2.  xhe  electric  field  is  enhanced  outside  the  auroral  are 
and  considerably  reduced  inside. 
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ELECTRIC  FIELD 


BIRKELAND  CURRENT 


DISTANCE  (Kml 

Fig.  5.4  — Electric  field,  Birkeland  current  density  and  anomalous  resistance 
of  1 cm^  flux  tube  for  the  current  driven  anomalous  resistance  mechanism. 
The  resistance  is  uniform  at  the  Spitzer  value  except  where  the  current  ex- 
ceeds the  instability  threshold  of  4.8  X 10~®  amp  m“2^  in  which  case  the 
resistance  increases  rapidly.  The  Birkeland  current  density  has  a square 
wave  character  with  the  amplitude  of  the  instability  current  threshold.  The 
electric  field  is  quite  similar  to  that  in  Fig.  5.3. 


has  the  shape  of  a square  wave  to  either  side  of  the  center  of  the  arc. 
The  resistance  of  a 1 cm'  flux  tube  is  also  shown.  We  see  that  it 


r 


peaks  at  the  steep  density  gradients  on  either  side  of  the  arc,  and  has 
a maximum  value  which  is  about  500  times  larger  than  the  Spitzer  resis- 
tivity . 

Figure  5-5  shows  the  results  of  a preliminary  attempt  to 
model  the  beam  driven  parametric  resistance  mechanism.  Papadopoulos 
and  Coffey®®  showed  that  the  resistance  is  roughly  proportional  to  the 
energy  in  the  precipitating  beam  of  electrons.  Here  we  assume, 
similarly,  that  the  resistance  is  proportional  to  the  total  ionization 
production  rate  in  the  auroral  arc.  This  leads  to  a resistance  which 
is  large,  0.5  esu,  at  the  center  of  the  arc  but  has  the  Spitzer  value 
immediately  outside.  The  results  show  that  both  the  Birkeland  current 
and  electric  field  peak  very  sharply  in  the  small  resistance  region 
immediately  outside  the  arc.  To  the  author's  knowledge,  such  strong 
peaking  has  not  been  observed  in  the  auroral  data. 

In  the  preceding  results  both  precipitating  and  cold  iono- 
spheric particle  currents  where  treated  identically.  In  all  these 
cases  the  Birkeland  current  was  balanced  across  the  center  of  the 
auroral  arc  density  enhancement  and  there  was  no  evidence  of  an  ener- 
getic precipitating  current.  The  next  three  samples  will  show  results 
in  which  the  effect  of  the  precipitating  particles  is  clearly  seen. 

In  these  cases  it  is  assumed  that  only  the  cold,  ionospheric  particle 
return  currents  experience  anomalous  resistance.  The  precipitating 
particles  carry  current  resistance-free. 
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Fig.  5.5  — Electric  field  and  Birkeland  current  density  for  the  beam  driven, 
parametric  resistance  mechanism.  Here  the  resistance  has  the  Spitzer  value 
outside  the  auroral  arc,  but  is  very  large  and  proportional  to  the  precipitat- 
ing energy  inside.  The  Birkeland  current  density  is  small  inside  the  arc  but 
becomes  large  at  the  edges.  The  electric  field  is  also  very  large  at  the  border 
of  the  arc. 
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Figure  5*6  shovjs  one  of  these  results  in  which  the  resistance 
is  large,  7 esu,  but  localized  to  the  auroral  arc  field  lines.  In  this 
case  the  precipitating  particle  current  (short  dash  line)  carries  almost 
all  of  the  Birkeland  current  inside  the  arc.  ITie  ionospheric  electron 
return  current  exists  outside  the  arc  in  a region  of  low  resistance. 
These  results  are,  likely,  a better  model  for  the  beam  driven  parametric 
resistance  mechanism  than  those  of  Figure  5-5- 

Figure  5*7  shows  a similar  result  but  in  this  case  the  resist- 
ance is  smaller,  O.J  esu,  and  is  constant.  Here,  the  precipitating 
current  is  shown  by  the  short  dash  line,  the  cold,  ionospheric  return 
current,  by  the  long  dash  line,  and  the  net  Birkeland  current,  by  the 
solid  line.  Again,  the  effect  of  the  precipitating  particles  is  clearly 
evident  in  the  total  Birkeland  current  density.  Nevertheless,  the 
ionospheric  return  current  is  substantial  inside  the  auroral  arc  and 
reduces  the  current  density  in  that  region.  This  effect  becomes  even 
more  apparent  in  the  results  of  Figure  5 *8 • For  this  result,  the 
constant  resistance  was  reduced  a factor  of  5 from  that  in  Figure  5-7- 
The  net  Birkeland  current  density  in  Figure  5-8  no  longer  shows  an 
extrema  in  the  region  of  the  arc  precipitation  and  is,  here  again, 
balanced  across  the  center  of  the  auroral  arc.  Any  further  reduction 
in  the  anomalous  resistance  would  lead  to  almost  complete  cancellation 
of  the  precipitating  particle  current  by  the  cold,  ionospheric  return 
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ELECTRIC  FIELD 


Fig.  5.6  — Electric  field  mid  Birkeland  current  density  for  resistance  free 
precipitating  particle  current.  The  resistance  to  the  ionospheric  return  cur- 
rent has  the  Spitzer  value  outside  the  auroral  arc  and  is  very  large  inside. 
The  short-dash  curve  shows  the  precipitating  current;  the  long-dash  curve 
shows  the  ionospheric  return  current;  and  the  solid  line  is  the  net  current 
density.  The  pre'cipitating  energetic  particles  carry  almost  all  of  the  electron 
current  into  the  ionosphere.  The  return  currents  are  strong  in  the  Spitzer 
resistance  region  outside  the  auroral  arc.  This  is  likely  a better  result  for 
the  beam  driven  resistance  mechanism  than  shown  in  Fig.  5.5. 
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ELECTRIC  FIELD 


Fig.  5.7  — Electric  field  and  Birkeland  current  density  for  a similar  case  to 
Fig.  5.6.  Here  the  resistance  is  spatially  uniform  and  has  a moderate  value, 
0.3  esu.  The  presence  of  the  resistance  free  precipitating  current  (short-dash 
curve)  is  evident  in  the  net  Birkeland  current  density  (solid  curve),  however, 
it  does  not  account  tor  more  than  about  half  of  the  net  downward  electron 
flow.  The  electric  field  is  similar  to  that  in  Figs.  5.3  and  5.4. 
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Fig.  5.8  — Results  for  the  same  conditions  as  in  Fig.  5.7  except  that  the  re- 
sistance is  reduced  a factor  of  5.  The  contribution  of  the  precipitating  cur- 
rent to  the  net  Birkeland  current  density  is  now  much  less  obvious.  The 
spatial  distribution  of  the  Birkeland  current  is  similar  to  that  in  the  first  3 
figures.  Further  reduction  of  the  anomalous  resistance  would  lead  to  a re- 
sult like  Fig.  5.1. 
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DISCUSSION 

In  the  preceding  section  a large  number  of  examples  were  pre- 
sented of  the  possible  effects  of  anomalous  resistivity  by  an  auroral 
Birkeland  current*  The  problem  which  remains  is  to  sort  out  the 
effects  and  to  develop  some  idea  of  their  possible  importance  and 
relevance.  That,  we  will  now  attempt  to  do. 

We  think  that,  taken  as  a whole,  the  results  presented  here 
are  not  in  great  disagreement  with  the  available  data  on  auroral 
Birkeland  currents  and  electric  fields.  That  is  to  say,  using  a model 
similar  to  that  presented  here  \drich  makes  reasonable  assumptions  con- 
cerning the  magnetospheric  electric  field  and  the  anomalous  resistivity 


mechanism  can  lead  to  a good  fit  between  the  model  results  and  the 


auroral  measurements.  We  therefore  feel  that  there  is  some  validity 

in  applying  the  model  results  to  the  auroral  currents  and  resistivity.  j 

i 

The  model  results  indicate  that  the  anomalous  resistivity  on  { 


auroral  field  lines  is  a factor  of  10^  to  2 x 19^  greater  than  the 
Spitzer  value.  There  are  a number  of  reasons  for  this.  If  the 
resistivity  is  too  small,  the  Birkeland  current  density  can  easily 
exceed  the  observed  range  of  the  measurements,  1 - 20  x 10”®  amp  m”'  . 
Moreover,  as  was  seen  in  the  results  of  Figure  5*3*  a small  resistance 
allows  the  ionospheric  return  current  to  almost  completely  neutralize 
the  precipitating  electron  current.  The  upper  value  of  the  resistivity 
is  set  by  plasma  instability  considerations  and  has  been  discussed 
elsewhere 

In  addition,  the  model  results  would  favor  a resistivity 
which  varies  slowly  in  space.  Rapid  changes  in  the  resistance  of 
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adjacent  flux  tubes  leads  to  sharp  extrema  in  both  the  Birkeland  cur- 
rents and  in  the  electric  field.  Spatial  changes  in  the  resistance  are 
important,  however,  as  they  greatly  influence  the  spatial  distribution 
of  the  Birkeland  current  density.  This  is  particularly  well  shown  in 
the  contrasting  results  of  Figures  “'.4  and  5 •5-  The  effect  of  a 
spatially  non-constant  resistance  is  to  alter  the  electric  field  so 
that  the  Birkeland  current  is  stronger  in  regions  of  lower  resistivity 
and  weaker  in  regions  of  higher  resistivity. 

Finally,  we  would  make  an  observation  regarding  the  accelera- 
tion of  auroral  particles  by  geomagnetic  field-aligned  potential  drops 
which  are  created  by  anomalous  resistivity.  For  the  results  sho\.m  here, 
the  largest  field-aligned  potential  drop  is  110  volts.  Even  if  the 
magnetospheric  electric  field  were  10  times  stronger,  it  could  not 
explain  the  energization  of  auroral  electrons  to  more  than  1 kilovolt. 
This  does  not  mean  that  the  mechanism  does  not  occur,  but  only  that  it 
cannot  occur  on  the  scale  size  (~  10  km)  of  a single  auroral  arc.  In 
fact  we  will  see  an  indication  that  it  does  occur  in  a more  realistic 
model  of  auroral  arcs  which  follows. 

5 .4  RESULTS  FOR  A MORE  REALISTIC  MODEL 

The  results  shown  and  discussed  in  the  previous  sections  were 
derived  from  a very  simple  model  of  an  auroral  arc  and  therefore  would 
not  be  expected  to  describe  accurately  the  fields  and  currents  in  an 
actual  aurora.  The  simple  model  was  used  in  order  to  show  more  clearly 
the  effects  of  anomalous  resistance  on  auroral  Birkeland  currents. 

Here  we  will  attempt  to  show  that  anomalous  resistivity  is  important  in 
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the  auroral  oval  by  use  of  a more  realistic  model.  ITris  model  was 
devised  in  an  attempt  to  understand  and  explain  the  OCO-6  electric 
field  measurements  from  the  midnight  region  of  the  auroral  oval  pre- 
sented by  Maynard. 

The  model  requires  the  assumption  of  the  magnetospheric 
electric  field  and  energetic  particle  precipitation  pattern  along  a 
north-south  meridian  through  the  premidnight  auroral  oval.  The  north- 
ward electric  field  chosen  is  l6  millivolts/meter,  constant  across  the 
oval;  and  no  northward  field  outside  of  the  oval.  The  westward 
electric  field  is  taken  to  be  25  millivolts/meter  north  of  the  oval, 
decreasing  linearly  to  zero  from  the  northern  to  southern  border  of 
the  oval,  and  zero  in  the  region  south  of  the  auroral  oval.  This 
choice  for  the  magnetospheric  electric  field  is  somewhat  unconventional 
and  controversial,  however  it  is  in  agreement,  we  think,  with  the  flow 
pattern  proposed  for  the  plasma  sheet  by  Rostoker  and  Bostrom®®  for 
the  immediate  premidnight  region. 

The  precipitation  pattern  chosen,  is  shown  as  the  short  dash 
line  in  Figure  -9 ■ It  features  a broad  (400  km  wide}  precipitation 
region  in  the  southern  and  central  regions  of  the  oval  and  a structured 
region  (containing  2 arcs  or  bands ) at  the  northern  border  of  the 
auroral  oval.  We  would  interpret  Lne  broad  region  as  plasma  sheet 
particle  precipitation  and  the  structured  region  as  inverted  - V 
precipitation.  The  precipitation  pattern  has  been  designed  to  be 
similar  to  the  Ariel  4 particle  data  of  Craven  and  Frank. 

Figure  5-9  also  presents  the  resulting  field  aligned  current 
systems  when  a uniform,  large  (7  esu^  anomalous  resistance  is  assumed. 


The  solid  line  shows  the  net  field-aligned  current  density,  the  long 
dash  line  shows  the  cold  ionospheric  electron  current,  and  the  short 
dash  line  shows  the  precipitating  current.  The  direction  and  strength 
of  the  net  current  is  consistent  with  the  TRIAD  satellite  date.^'^^  The 
important  effect  of  the  anomalous  resistance  is  to  cause  the  hot  and 
cold  particle  currents  to  oppose  each  other  in  the  broad  precipitation 
region  and  to  have  the  same  direction  in  the  structured  precipitation 
region.  Rocket  magnetometer  data  have  been  interpreted  to  indicate 
that  both  hot  and  cold  electrons  flow  downward  into  arcs  T see  Anderson 
and  Cloutier^®"^  and  references  therein]. 

In  addition  to  the  cold  ionospheric  current  the  long  dash 
line  in  Figure  5*9  is  a plot  of  the  field-aligned  potential  difference. 
The  units  in  this  case  are  kilovolts  . In  the  structured  region  the 
potential  is  directed  upward  and  becomes  as  large  as  2 kilovolts. 
Presumably  this  potential  could  serve  to  accelerate  electrons  into  the 
ionosphere.  In  the  broad  precipitation  region  the  potential  is  doti/nward 
and  thus  would  discourage  electron  precipitation  and  encourage  ion 
precipitation . 

Figure  5*10  is  a comparison  between  a measured  electric 
field and  the  model  electric  field  result.  The  antenna  for  the 
measurement  was  oriented  north-east  south-west  and  therefore  the  model 
field  was  plotted  on  the  same  axis  . In  each  case  the  major  features  of 
the  field  are  the  same.  At  the  southern  border  the  field  is  northward; 
it  passes  through  zero  at  the  Harang  Discontinuity  (Maynard's  identifi- 
cation' and  becomes  northward  through  the  central  region  of  the  auroral 
oval.  At  the  northern  border,  both  fields  oscillate  rapidly  before 


r 
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turning  basically  westward  in  the  polar  region,  'llie  agreement  between 
the  model  electric  field  and  the  measurement  is  good. 

In  this  paper  we  have  shown  that  a model  incorporating 
anomalous  resistivity  to  field  aligned  currents  can  provide  a descrip- 
tion of  auroral  electric  fields  and  currents,  which  is  in  reasonable 
agreement  with  auroral  measurements.  It  was  also  shown  that  the 
effects  of  the  resistivity  were  important  to  the  Birkeland  current 
densities  derived  therefrom.  A better  evaluation  of  resistivity  and  its 
effect  on  auroral  fields  and  currents  can  be  made  only  by  incorporating 
the  actual  plasma  mechanisms  and  their  nonlinear  behavior  in  the  auroral 
model. 


I 


I 

5-24 


L 


FIELD-ALIGNED  CURRENT 


SOUTH  DISTANCE  (km)  NORTH 

Fig.  5.9  — The  field-aligned  current  density  for  a more  realistic  model  of  auroral 
arcs  in  the  premidnight  auroral  oval.  The  short-dash  curve  shows  the  precipitating 
energetic  electron  current.  It  is  composed  of  a broad  precipitation  region  (plas- 
ma sheet  particles)  in  the  southern  and  central  part  of  the  oval,  and  two  discrete, 
narrow  precipitation  structures  (auroral  arcs  or  hands)  at  the  northern  Irorder  of 
the  oval.  The  long-dash  curve  shows  the  cold-ionospheric  electron  current.  In 
the  southern  and  central  portions  of  the  oval  the  hot  and  cold  particle  currents 
oppose  each  other,  while  at  the  northern  border  they  are  m the  same  direction. 
Tiie  solid  curve  shows  the  net  field-aligned  current  density  which  results.  In 
addition  to  the  cold  ionospheric  current,  the  long-dash  curve  shows  the  potential 
drop  along  the  geomagnetic  field;  the  appropriate  units  are  kilovolts.  The  poten- 
tial is  upward  and  accelerates  precipitating  electrons  into  the  discrete  precipitation 
region,  while  opposing  electron  precipitation  in  the  broad  region. 
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Fig.  5.10—  A comparison  between  a measured  electric  field  [Maynard, 
1974]  and  the  model  electric  field  for  the  currents  of  Fig.  5.9.  Both 
the  measured  and  model  fields  have  the  same  features.  The  field  is 
northward  at  the  southern  border  of  the  auroral  oval,  becomes  south- 
westward  through  the  central  oval,  and  sho’-s  similar  structure  in  the 
discrete  precipitation  region  at  the  northern  border. 


Section  6 
SUMMARY 

In  the  introduction,  we  have  pointed  out  the  considerable 
similarities  between  the  effects  of  a HANE  and  those  of  an  auroral 
disturbance.  These  similarities  include  precipitation  of  energetic 
particles,  generation  of  field-aligned  current  systems,  atmospheric 
heating  and  generation  of  winds,  radiation  throughout  the  electro- 
magnetic spectrum,  along  with  other  plasma  phenomena.  One  result  of 
these  various  processes  is  the  creation  of  communications  disturbances 
and  interruptions  which  can  be  far  reaching  and  of  long  duration.  It 
is  the  study  of  these  phenomena  and  their  resulting  effect  on  communica- 
tions systems  which  is  the  motivation  for  the  work  renorted  here. 

This  is  not  to  say  that  the  aurora  demonstrates  these  nuclear 
phenomena  in  the  same  magnitude  and  to  the  same  extent  as  a HANE,  but 
rather  that  the  same  physical  laws  and,  therefore,  similar  models 
describe  the  effects  in  both  events.  The  aurora  can  therefore  be  used, 
in  the  absence  of  nuclear  testing,  as  a benchmark  for  nuclear  phenom- 
enology. Moreover,  the  physical  models  developed  in  the  auroral 
program  can  be  applied,  with  some  modification,  to  the  study  of  nuclear 
phenomenology.  The  auroral  program,  in  this  manner,  is  an  integral 
part  of  the  study  of  HANE  effects  on  communication  systems. 

A program  for  developing  an  auroral  model  must  address  three 
fundamental  technical  areas:  1)  the  deposition  of  energy  in  the  upper 
atmosphere  by  energetic  electrons  and  protons  2)  the  chemical  reactions 
initiated  by  the  deposited  energy  and  5)  the  electrodynamical  aspects 
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such  as  the  structures  of  fields  and  currents.  This  program  closely 
parallels  the  approach  necessary  in  modeling  HANE  effects  and  provides 
valuable  insights  for  the  understanding  of  HANE  phenomenology  as  well 
as  serving  as  a measuring  standard  for  our  ability  to  model  perturbations 
(or  HANEs)  in  the  upper  atmosphere. 

During  the  past  year,  significant  improvements  have  been  in- 
corporated into  the  NRL  proton  deposition  model.  The  number  of  atmo- 
spheric constituents  in  the  model  has  been  increased  with  the  addition 
of  0 and  0^  to  the  N already  in  the  program.  The  photon  energy 
degradation  treatment,  based  on  the  continuous  slowing  down  approxima- 
tion, has  been  improved  by  the  inclusion  of  the  influence  of  the 
geomagnetic  field  on  the  particle  energy  loss. 

Progress  has  also  been  made  in  electron  energy  deposition. 

The  NRL  electron  transport  code  has  predicted  higher  backscatter  than 
other  models  because  of  the  way  the  energy  loss  process  is  treated.  A 
detailed  study  reveals  high  sensitivity  of  the  backscattered  flux  to 
variations  in  the  large  angle  scattering  cross  sections.  The  numerical 
solutions  to  the  equation  of  transfer  are  shown  to  be  accurate  to  at 
least  three  significant  figures. 

Auroral  chemistry  studies  were  focused  on  the  large  enhance- 
ments of  NO  concentration  in  auroral  arcs  inferred  from  observations. 

An  analysis  in  terms  of  a model  incorporating  the  possible  channels  of 
NO  buildup  has  indicated  that  NO  enhancements  in  excess  of  10^  of  the 

concentration  of  0 would  be  improbable. 

2 

Improvements  have  been  made  in  the  auroral  electrodynamics 
models.  A detailed  assessment  of  the  importance  of  anomalous  plasma 


6-2 


resistivity  in  modeling  auroral  electric  fields  and  currents  has  been 
made.  These  model  results  have  been  shown  to  be  in  good  agreement  with 
auroral  measurements. 

The  studies  above  have  all  contributed  to  other  areas  of  the 
nuclear  effects  program.  For  instance,  the  proton  energy  deposition 
method  has  been  modified  to  describe  the  deposition  of  nuclear  debris 
ions  and  has  given  results  which  compare  well  with  data  from  bomb  test 
programs.  It  is  now  used  as  an  important  part  of  HANE  simulations. 
Changes  have  also  been  made  in  HANE  chemistry  codes  as  a direct  result 
of  auroral  chemistry  studies.  The  auroral  electrodynamics  studies  and 
models  show  promise  of  making  Important  contributions  in  the  near 
future  as  nuclear  phenomenology  studies  are  carried  to  later  times. 

Much  of  the  theoretical  developments  in  the  auroral  models  are  directly 
applicable  to  the  electrostatic  evolution  of  the  nuclear  debris  and 
ionized  plasma  which  is  of  great  importance  in  evaluating  very  late 
time  effects  of  HANEs  on  defense  satellite  communications  systems. 

Future  work  will  be  devoted  to  studying  the  occurrence  and 
development  of  plasma  irregularities  in  the  auroral  ionosphere.  The 
electrodynamic  models  used  for  this  study  are  closely  related  to  the 
late-time  electrostatic  models  of  HANEs.  It  is  expected  that  these 
studies  will  allow  a direct  connection  to  be  made  between  scintilla- 
tions at  auroral  latitudes  observed  in  the  Wideband  Satellite  Experi- 
ment and  those  which  would  occur  in  the  nuclear  environment.  The 
auroral  modeling  program  is  unique  in  allowing  a direct  connection  to 
be  made  between  naturally  occurring  scintillations  and  HANE  effects  on 


communications . 


In  conclusion,  we  point  out  that  the  auroral  modeling  program 
has  made  direct,  important  contributions  to  llANE  phenomenology  studies. 
Moreover,  it  is  expected  that  such  benefits  will  continue  to  accrue. 

As  interest  increases  in  nuclear  disturbance  effects  at  later  times, 
the  auroral  modeling  effort  and  the  electrostatic  methods  and  techniques 
developed  through  it  will  become  a foundation  of  future  developments  in 
HANE  phenomenology.  In  addition  it  will  allow  a connection  to  be  made 
between  nuclear  weapons  effects  on  communications  and  measurements 
which  are  conducted  in  the  natural  environment. 
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OICY  ATTN  CODE  WA501  NAVY  NUC  PRGMS  OFF 


DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON,  D.C.  20376 
OICY  ATTN  NSP-2141 

OICY  ATTN  NSSP-2722  FRED  WIMBERLY 


NAVAL  SPACE  SYSTEM  ACTIVITY 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALIF.  90009 

OICY  ATTN  A.  B.  HAZ2ARD 


COMMANDER 

ADC/DC 

ENT  AFB,  CO  80912 

OICY  ATTN  DC  MR.  LONG 


COMMA, NOE  R 
ADCOM/XPD 
ENT  AFB.  CO  80912 

OICY  ATTN  XPQDQ 


AF  GEOPHYSICS  LABORATORY,  AFSC 
HANSCOM  AFB,  MA  01751 

OICY  ATTN  OPR  HAROLD  GARDNER 

OiCY  ATTN  OPR  JAMES  C.  ULWICK 

OICY  ATTN  LKB  KENNETH  5.  W.  CHAMPION 

OICY  ATTN  OPR  ALVA  T.  STAIR 

OICY  ATTN  SUOL  RSCH  LIB 

OICY  ATTN  PHP  JULES  AARONS 

OICY  ATTN  PHD  JURGEN  BUCHAU 

OICY  ATTN  PHD  JOHN  MULLEN 


AF  WEAPONS  LABORATORY,  AFSC 
KiRTLAND  AFB,  NM  87117 


0 ICY 

ATTN 

SUL 

0 ICY 

ATTN 

CA 

ARTHUR  H.  GUENTHER 

0 ICY 

ATTN 

DYC 

CAPT  L.  WITTWER 

OICY 

ATTN 

SAS 

JOHN  M.  KAMM 

0 ICY 

ATTN 

DYC 

CAPT  MARK  A.  FRY 

AFTAC 

PATRICK  AFB,  FL  32925 

OICY  ATTN  TF/MAJ  WILEY 
OIC^  ATTN  TN 


AIR  FORCE  AVIONICS  LABORATORY,  AFSC 
WR I GHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AFAL  AAB  H.  M.  HARTMAN 
OICY  ATTN  AAD  ALLEN  JOHNSON 


HEADQUARTERS 

electronic  systems  DIVISION/XR 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  XRC  LTC  J.  MORIN 
OICY  ATTN  XRE  LT  MICHAELS 


HEADQUARTERS 

electronic  systems  DIVISION/YS 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  YSEV 
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f 


M 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WR I GHT-PATTERSON  AF3,  OH  45A33 
OICY  ATTN  NICD  LIBRARY 
OICY  ATTN  ETD  B.  L.  BALLARD 


HQ  USAF/RD 

WASHINGTON,  D.C.  20330 
OICY  ATTN  RDQ 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13440 

OICY  ATTN  EMTLD  DOC  LIBRARY 
OICY  ATTN  OCSE  V.  COYNE 


SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

OICY  ATTN  SZJ  MAJOR  LAWRENCE  DOAN 


COMMANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB,  NB  68113 

OICY  ATTN  XPFS  MAJ  BRIAN  G.  STEPHAN 
OICY  ATTN  ADWATE  CAPT  BRUCE  BAUER 
OICY  ATTN  NRT 


HEADQUARTERS 

electronic  systems  division  (AFSC) 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  JIM  DEAS 


SAMSO/ YA 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

OICY  ATTN  YAT  CAPT  L.  BLACKWELDER 


SAMSO/SK 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 


OICY  ATTN  SKA  LT  MARIA  A.  CLAVIN 


SAMSO/MN 

NORTON  AFB,  OA  92419 
(MINUTEMAN) 

OICY  ATTN  MNNL  LTC  KENNEDY 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01733 

OICY  ATTN  ETEI  A.  LORENTZEN 


U.5.  ENERGY  RSCH  AND  DEV  ADMIN 


f 


EG6G,  INC. 

LOS  alamos  division 

P.  0.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  JAMES  R.  BREEDLOVE 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  TECH  INFO  DEPT  L-3 
OICY  ATTN  RONALD  L.  OTT  L-531 
OICY  ATTN  DONALD  R.  DUNN  L-156 
OICY  ATTN  RA5PH  S.  HAGER  L-31 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1663 

LOS  alamos,  NM  87545 

OICY  ATTN  DOC  CON  FOR  ERIC  LINDMAN 

OICY  ATTN  DOC  CON  FOR  R.  F.  TASCHEK 

OICY  ATTN  DOC  CON  FOR  ERIC  JONES 

OICY  ATTN  DOC  LON  FOR  JOHN  S.  MAL I< 

OICY  ATTN  DOC  CON  FOR  MARTIN  TIERNEY  J-IO 

OICY  ATTN  DOC  CON  FOR  JOHN  2INN 


SANDIA  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  J.  P.  MARTIN  ORG  1732 

OICY  ATTN  DOC  CON  FOR  W.  D.  BROWN  ORG  1353 

OICY  ATTN  DOC  CON  FOR  A.  DEAN  THORN3ROUGH  ORG  1245 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  A.  DAHLGREN  ORG  1722 


OTHE*?  government 


DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TELECOMMUNICATIONS 
INSTITUTE  FOR  TELECOM  SCIENCE 
BOULDER,  CO  80302 

OICY  ATTN  WILLIAM  F.  UTLAUT 
OICY  ATTN  G.  REED 


NATIONAL  OCEANIC  & ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  JOSEPH  H.  POPE 
OICY  ATTN  RICHARD  GRUBB 
OICY  ATTN  C.  L.  RUFFNACH 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 

OICY  ATTN  ATS-6  OFC  P.  CORRIGAN 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 


P.  0.  BOX 

929  57 

LOS  ANGELES,  CA 

90009 

OICY 

ATTN 

IRVING 

M.  GARFUNKEL 

OICY 

ATTN 

T 

. M . 

SALMI 

0 ICY 

ATTN 

V 

. JOSEPHSON 

0 ICY 

ATTN 

S 

. P . 

BOWER 

0 ICY 

ATTN 

N 

. D . 

5TOCKWELL 

OICY 

ATTN 

P 

. P . 

OLSEN  120  RM  2224E 

0 ICY 

ATTN 

F 

. E . 

BOND  RM  5003 

0 ICY 

ATTN 

J 

. E . 

CARTER  120  RM  2209 

OICY 

ATTN 

F 

. A . 

MORSE  A6  RM  2407 

analytical  SYSTEMS  ENGINEERING  CORP 

5 OLD  CONCORD  ROAD 

BURLINGTON,  MA 

01803 

OICY  ATTN 

RADIO  SCIENCES 

BOEING  COMPANY, 

THE 

P.  0.  BOX  3707 

SEATTLE,  WA  98124 

OICY  ATTN 

GLEN  KEISTER 

OICY  ATTN 

D.  MURRAY 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
3175  MIRAMAR  ROAD 
LA  JOLLA,  CA  92057 

OICY  ATTN  HENRY  G.  BOO<ER 


BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  research  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMEO  A.  DELI8ERIS 


CHARLES  STARK  DRAPER  LABORATORY,  INC 
555  technology  square 
CAMBRIDGE,  MA  02159 

OICY  ATTN  D.  B.  COX 

OICY  ATTN  J.  P.  GILMORE  MS  63 
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COMPUTER  SCIENCES  CORPORATION 
P.  0.  BOX  530 
6565  ARLINGTON  BLVD 
falls  church,  VA  22046 
OICY  ATTN  H.  BLANK 
OICY  ATTN  JOHN  SPOOR 


COMSAT  laboratories 


LINTHICUM 

ROAD 

CLARKSBURG 

, MD 

20  7 34 

OICY 

ATTN 

R.  R.  TAUR 

CORNELL 

UNI 

VERS  ITY 

DEPARTMENT 

OF  electrical  ENGI 

NEERING 

I THACA, 

NY 

148  50 

0 ICY 

ATTN  D.  T.  FARLEY 

JR 

ESL 

INC  . 

495 

JAVA 

DRIVE 

SUNNYVALE 

, CA 

94086 

OICY 

ATTN 

R.  K. 

STEVENS 

OICY 

ATTN 

J.  ROBERTS 

OICY 

ATTN 

JAMES 

MARSHALL 

OICY 

ATTN 

V . L . 

MOWER 

0 ICY 

ATTN 

C . W . 

PRETTIE 

FORD  AEROSPACE  6 COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  alto,  CA  94303 

OICY  ATTN  J.  T.  MATTINGLEY  MS  X22 


GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  80RTNER  SPACE  SCI  LAB 
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GENE’^Al  electric  company 

TEMPO-CENTER  FOR  ADVANCED  STUDIES 
81S  STATE  STREET  CP.O.  DRAWER  OQ) 
SANTA  BARBARA,  CA  93102 
OICY  ATTN  OASIAC 
OICY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHENS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  william  MCNAMERA 
OICY  'ATTN  B.  GAMBILL 
OICY  ATTN  MACK  STANTON 


general  RESEARCH  CORPORATION 

P.  0.  BOX  3587 

SANTA  BARBARA,  CA  95105 

OICY  ATTN  JOHN  ISE  JR 
OICY  ATTN  JOEL  GAR3ARIN0 


geophysical  institute 

UNIVERSITY  OF  ALASKA 

FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

OICY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

OICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 


GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  OIV 
77  A STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  J.  M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFHARD 
OICY  ATTN  JOEL  8ENGSTEN 


HARRIS  CORP,  E.S.O. 

P.  0.  BOX  37 
MELBOURNE,  FL  32901 

OICY  ATTN  ADV  PROG  DEPT  OR.  CARL  DAVIS 


HSS,  [NC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01750 

OlCr  ATTN  DONALD  HANSEN 


INTL  TEL  S TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NU  07110 

OICY  ATTN  technical  LIBRARY 


JAY  COR 

1401  camino  del  mar 

DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
applied  physics  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MO  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 


LOCKHEED  MISSILES  & SPACE  CO  INC 
P.  0.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 


LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 

OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 

OICY  ATTN  BILLY  M.  MCCORMAC  DEPT  52-54 


KAMAN  SCIENCES  CORP 

P.  0.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 

OICY  ATTN  B.  J.  BITTNER 


LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 
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M.l.T.  LINCOLN  LABORATORY 
P.  0.  BOX  75 
LEXINGTON,  MA  02173 

OICY  ATTN  LIB  A-082  FOR  DAVID  M.  TOWLE 
OlCY  ATTN  MR.  WALDEN  X115 
OICY  ATTN  JAMES  H.  PANNELL  L-246 
OICY  ATTN  D.  CLARK 


MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  BILL  OLSON 


MISSION  RESEARCH 
735  STATE  STREET 


CORPORAT I ON 


BARBARA,  i 

CA 

93101 

0 ICY 

ATTN 

P . 

F I SCHER 

0 ICY 

ATTN 

W . 

F.  CREVIER 

0 ICY 

ATTN 

STEVEN  L.  GUTSCHE 

0 ICY 

ATTN 

0. 

SAPPENF lELD 

0 ICY 

ATTN 

R . 

BOGUSCH 

0 ICY 

ATTN 

R . 

HENDR I CK 

0 ICY 

ATTN 

RALPH  KILB 

OICY 

ATTN 

DAVE  SOWLE 

OICY 

ATTN 

F . 

FAJEN 

OICY 

ATTN 

M . 

SCHEI BE 

0 ICY 

ATTN 

CONRAD  L.  LONGMIRE 

OICY 

ATTN 

WARREN  A.  SCHLUETE 

MITRE  CORPORATION, 
P.  0.  BOX  208 
BEDFORD,  MA  01730 


THE 


OICY  ATTN  J.  C.  KEENAN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  CHIEF  SCIENTIST  W, 
OICY  ATTN  S.  A MORIN 
OICY  ATTN  C.  E.  CALLAHAN 


SEN 


PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 

0 ICY  ATTN  F . C . FI  ELD  UR 


I 
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PHOTOM.ETRICS,  INC. 

442  MAriRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 


P^^YSICAL  dynamics  INC. 

P.  0.  SOX  21589 
SEATTLE,  WA  98111 

OICY  ATTN  E.  J.  FREMOUW 


PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  94701 

OICY  ATTN  A.  THOMPSON 
OICY  ATTN  JOSEPH  B.  WORKMAN 


R 6 0 ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  RICHARD  LATTER 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 


SCIENCE  applications,  INC. 

P . 0 . BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGAL 


CO. 

5 28  BOSTON  ="OST  ROAD 
S.'DB-R'',  ma  01776 

01  Y A’-TN  BARBARA  ADAMS 


SC!£\CE  APPl  I CA  MO.'.S  , I ^JC  . 
HUN^SV  1 ^LE  D W I S 1 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE.  Al  55805 

OICY  ATTN  dale  H.  DIVIS 


SCIENCE  APPL I CAT ! ON5,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

0 ICY  ATTN  B . ADAMS 


STANFORD  RESEARCH  INSTITUTE 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 

OICY  ATTN  ALAN  BURNS 

OICY  ATTN  G.  SMITH 

OICY  ATTN  L.  L.  COBB 

OICY  ATTN  DAVID  A.  uOHNSON 

OICY  ATTN  WALTER  G.  CHESNUT 

OICY  ATTN  CHARLES  L.  RING 

OICY  ATTN  WALTER  JAYE 

OICY  ATTN  M.  BARON 

OICY  ATTN  RAY  L.  LEADABRAND 


SYSTEM  DEVELOPMENT  CORPORATION 
4130  LINDEN  AVENUE 
DAYTON,  OH  45432 

OICY  ATTN  F.  G.  MEYER 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.  P.  BOQUIST 


TRW  DEFENSE  6 SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLE8UCH  Rl-2078 
OICY  ATTN  ROBERT  M.  WEBB  Rl-1150 


VISIOYNE,  INC. 

19  third  avenue 

north  west  industrial  park 

BURLINGTON,  MA  01805 

OICY  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.  W.  CARPENTER 
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